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Abstract 

A quantum kinetic formalism is developed to study the dynamical interplay of quantum and 



> 

in 

statistical-kinetic properties of non-equilibrium multi-parton systems produced in high-energy QCD 
processes. The approach provides the means to follow the quantum dynamics in both space-time and 

in ■ 

. energy-momentum, starting from an arbitrary initial configuration of high-momentum quarks and 

gluons. Using a generalized functional integral representation and adopting the 'closed-time-path' 



Green function techniques, a self-consistent set of equations of motions is obtained: a Ginzburg- 
| Landau equation for a possible color background field, and Dyson-Schwinger equations for the 

2-point functions of the gluon and quark fields. By exploiting the 'two-scale nature' of light-cone 
dominated QCD processes, i.e. the separation between the quantum scale that specifies the range 
of short-distance quantum fluctuations, and the kinetic scale that characterizes the range of statis- 
tical binary interactions, the quantum- field equations of motion are converted into a corresponding 
set of 'renormalization equations' and 'transport equations'. The former describe renormalization 
and dissipation effects through the evolution of the spectral density of individual, dressed partons, 
whereas the latter determine the statistical occurrence of scattering processes among these dressed 
partons. The renormalization equations and the transport equations are coupled, and hence must be 
solved self-consistently. This amounts to evolving the multi-parton system, from a specified initial 
configuration, in time and full 7-dimensional phase-space, constrained by the Heisenberg uncertainty 
principle. This quantum-kinetic description provides a probabilistic interpretation and is therefore 
of important practical value for the solution of the dynamical equations of motion, suggesting for 
instance the possibility of simulating the multi-particle dynamics with Monte Carlo methods. 

klaus@suryall.cern.ch 

CERN-TH. 95-181, July 1995 



1. INTRODUCTION 



In this paper I attempt to formulate an approach towards a fundamental and consistent 
description of the statistical properties of non-equilibrium quantum systems produced in 
high-energy QCD processes, which allows to follow the quantum dynamics in time and 
complete phase space starting from an initial configuration. It provides a flexible framework 
for a systematic analysis of typical problems associated with the quantum dynamics of such 
systems, including, e.g., multi-particle transport phenomena of gluons, quarks and hadrons, 
or, critical dynamics of phase-transition phenomena and spontanous symmetry breaking, 
or, quantum dissipation, entropy generation and multi-particle production. 

More specifically, the intentions are aimed towards a practically applicable description 
of the space-time evolution of a general initial system of gluons and quarks, characterized 
by some large energy or momentum scale, that expands, diffuses and dissipates according to 
the self- and mutual- interactions, and eventually converts dynamically into excited hadronic 
matter and a final state hadron system by a "phase transition". This scenario frames a 
wide class of QCD processes of both fundamental and phenomenological interest. For 
instance, the evolution of parton showers in the mechanism of parton-hadron conversion 
in elementary high-energy processes (e + e~-annihilation into hadrons, deep-inelastic lepton- 
nucleon scattering, or non-diffr active hadronic collisions), or, the description of formation, 
evolution and freeze-out of a quark- gluon plasma in ultra-relativistic heavy- ion collisions, 
or, the study of the dynamics of the QCD phase-transition from the deconfined, high- 
temperature partonic phase to a low temperature hadronic phase with the simultanous 
breakdown of chiral symmetry and the condensation of gluons and quarks in the vacuum, 
as it occurred during the early evolution of the Universe. 

In the present paper I will confine myself to the first stage, the high-energy quark-gluon 
phase, and develop a quantum kinetic formalism that allows to describe both the dissi- 
pative and dispersive dynamics of a multi-parton system in real time. This description is 
exclusively based on the fundamental QCD Lagrangian and its firmly established princi- 
ples. The second stage, the parton-hadron conversion and phase transition, on the other 
hand, requires supplementary phenomenological input to model the details of the confine- 
ment mechanism that are not known at present jj]]. Such a phenomenological approach to 
the real-time dynamics of parton-hadron conversion that models the transition within an 
effective field theory description has been proposed recently in Ref. Q. It is preferable, 
however, to keep the fundamental description of the first stage distinct from the less under- 
stood phenomenological aspects of the second stage, and therefore I will address the latter 
in a separate paper. 
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In general, the study of a high-energy multi-particle system and its quantum dynamics 
involves three essential aspects: first, the aspect of space-time, geometry and the structure 
of the vacuum; second, the quantum field aspect of the particle excitations; and third, 
the statistical aspect of their interactions. These three elements are generally intercon- 
nected in a non-trivial way by their overall dynamical dependence. Therefore, in order 
to formulate a quantum description of the complex non-equilibrium dynamics, one needs 
to find a quantum-statistical and kinetic formulation of field theory that unifies the three 
aspects self-consistently. With this paper I take steps towards this goal by combining three 
corresponding theoretical methods, namely, first, the closed-time-path (CTP) formalism 
j|, [|, ||, ||, |?], H (for treating initial value problems of irreversible systems), second, the non- 
local source theory ||, |llj (for incorporating quantum fluctuations), and third, transport 
theory based on Wigner function techniques (l^, [l3|, 14] (for a kinetic description of inho- 
mogenous non-equilibrium systems). In principle, a dynamical theory of non-equilibrium 
multi-particle systems as the above mentioned, should be described by an exact quantum 
kinetic theory of QCD. Over the past 10 years, elaborate works |l5], [l6|, |l7]] have put great 
effort into deriving a general QCD transport theory rigorously from first principles. Un- 
fortunately, due to a number of unresolved problems arising from the complexities of the 
non-abelian gauge structure of QCD, the derived gauge-covariant formalism remains an 
academic theory up to date. It is of little practical value, unless it is boiled down to the 
quasi-classical limit by a series of approximations yielding a mean-field description, which 
however cannot describe the production of physical particles and their spectra. 

I am less ambitous here in what concerns the generality, and instead put emphasis on ap- 
plicability to realistic physical situations, in particular to the type of lightcone-dominated 
processes that I classified above. This class of high-energy processes allows a clear dis- 
tinction between a short-distance quantum field theoretical scale and a larger distance 
statistical-kinetic scale. When described in a reference frame, in which the particles move 
close to the speed of light, the effects of time dilation and Lorentz contraction separate the 
intrinsic quantum motion of the individual particles from the statistical correlations among 
them. On the one hand, the quantum dynamics is determined by the self-interactions of 
the bare quanta, and by the possible presence of a coherent background field (or mean 
field in the Hartree-Fock sense), in case one desires to go beyond a description in the pure 
vacuum. This requires a fully quantum theoretical analysis including renormalization. On 
the other hand, the kinetic dynamics can well be described statistical-mechanically by the 
motion of the quasi-particles which arise from the 'dressing' of the bare quanta by their 
self-interactions and by the background field, plus the binary interactions between these 
quasi-particles. Such a distinct description of quantum and kinetic dynamics is possible, 
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because the quantum fluctuations are highly concentrated around the light cone, occurring 
at very short distances, and decouple to very good approximation from the kinetic evolu- 
tion which is dictated by comparably large space-time scales. As mentioned, the natural 
two-scale separation is just the consequence of time dilation and Lorentz contraction, and 
is true for any lightcone dominated process. In fact, at asymptotic energies the quantum 
fluctuations are exactly localized on the lightcone, and so the decoupling becomes perfect. 
This observation is the key to formulate a quantum kinetic description in terms of par- 
ticle phase-space densities, involving a simultanous specification of momentum space and 
space-time, because at sufficiently high energy, the momentum scale Ap of the individual 
particles' quantum fluctutions and the scale Ar of space-time variations of the system of 
particles satisfiy ApAr » 1, consistent with the uncertainty principle. 

With this physical input and utilizing the aforementioned theoretical tools, the analysis 
proceeds as follows. In the first step, covered in Sec. 2, I obtain, starting from the QCD 
Lagrangian, the CTP generating functional for the gluon and quark Green functions, being 
defined on a closed-time contour and incorporating initial state correlations. From the as- 
sociated effective action, one gets the quantum dynamical equations of motion, which are 
the CTP version of the Ginzburg-Landau equation and the Dyson-Schwinger equations. In 
the second step, described in Sec. 3, I make the transition from quantum field description 
to kinetic theory, by exploiting the two-scale nature of lightcone dominance, and more- 
over, choosing a ghost-free axial gauge for the gluon fields. As a result one obtains from 
the Dyson-Schwinger equations a set of kinetic equations, consisting of a renormalization 
equation, that describes the quantum dynamics in terms of short-distance self-interactions 
of gluons and quarks, plus a transport equation that describes the kinetic dynamics of re- 
laxation and collision processes in terms of the statistical interactions of the renormalized, 
dynamically dressed partons among each other. The renormalization equation and the 
transport equations are coupled, and hence must be solved self-consistently. This amounts 
to evolving the system under consideration from its initial configuration simultanously in 
position- and momentum-space, constrained by the Heisenberg uncertainty principle. Fi- 
nally, Sec. 4 closes with some concluding remarks, and the Appendices summarize, for each 
of the above aspects, the technical details which are only indicated in the text. 

The main findings can be summarized as follows. The dynamics of high-energy multi- 
parton systems can, under reasonable conditions, be described in a semi-classical manner: 
the partons can be considered as dressed quanta with a dynamical substructure and a 
corresponding form factor arising from the self-interactions. The space-time evolution of a 
system of many such dressed partons is then governed by their propagation along classical 
trajectories and mutual binary collisions, as determined by their density, cross-sections and 
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by quantum statistics. This emerging picture is of great practical value for formulating a 
systematic calculation scheme — in a sense the space-time generalization of the 'jet calculus' 
fL8| , |i~9| ]. In Sec. 3.6, I outline such a scheme. One of the greatest advantages of this kinetic 
description is that it provides a probabilistic interpretation of the time evolution in full 
7-dimensional phase-space, which suggests the opportunity to simulate the multi-particle 
dynamics as sequential Markov processes with Monte Carlo methods. 

Finally let me comment on placing this work in relation to existing literature. 

(i) The general ideas and techniques of the CTP functional integral formalism have orig- 
inally been introduced mainly by Schwinger (3|, Keldysh jjj, and by Kadanoff and 
Baym [||], more than thirty years ago. The most extensive review that sums up 
the current state of the art is probably the work of Chou et al. ||, with diverse 
exemplification of the wide class of physics applications. Further pedagogically ex- 
cellent presentations have been published by Calzetta and Hu Wi, and by Rammer 
and Smith ||. In the particular field of relativistic nuclear physics, the concepts 
have been pragmatically applied , e.g., by Li and McLerran f20fl, and by Zhang and 



Wilets [21|. Important contributions of fundamental studies have been made in the 
last years by Danielewicz |22[], and Mrowczinsky and Heinz |2^]. The goal to establish 
a quantum kinetic theory for QCD was pioneered by the ambitious efforts of Elze, 
Gyulassy, Heinz, and Vasak p5l [l6|, |l7]], which resulted in a rigidly general, gauge 
covariant formalism. However, the prize to pay is an intractable complexity that, 
without specific physics input, is essentially of aestethic value without much practical 
use. The new achievement of the present work from this perspective may be stated 
as the adaption of the general CTP formalism, applied to QCD, but with focus on 
situations where the multi-parton dynamics is characterized by a large energy scale 
and can be described reliably within perturbation theory in a physical gauge. 

(ii) The most related recent works from the viewpoint of attempting to tackle evolution 
of multi-parton systems at high energy are probably the innovative works of McLer- 



ran, Venugopalan, et al. [24], and of Makhlin [25], in which the issue of calculating 
parton distributions in the context of ultra-relativistic nuclear collisions is addressed. 
The former authors use a classical non-abelian field description of QCD to compute 
coherent initial state properties of colliding large nuclei, whereas the latter focusses 
on a quantum field description of final state correlations of produced particles. Both 
approaches however do not attempt to address explicitly the space-time evolution of 
the multi-parton ensemble emerging from the nuclear collision, which is the main goal 
of the present paper. 
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(iii) The key elements to address the space-time evolution are provided by the Wigner func- 
tion techniques, which date back to Wigner's work on transport phenomena and 
are reviewed in in e.g. [13, 14]. Although widely exploited in condensed matter and 



plasma physics, these tools for quantum kinetics of many-body systems have hardly 
been applied to describe high energy non-equilibrium dynamics in QCD. New in the 
present work is the synthesis of quantum dynamics on the basis of the renormalization 
group of QCD, and quasi-particle kinetics within relativistic transport theory. The 
combination of these two aspects forms the foundation of the self-consistent treat- 
ment that entails a thorough consideration of the renormalization problem, which is 
commonly avoided in other applications (see however Ref. 0]). 

(iv) The machinery of perturbative QCD for light-cone dominated high energy processes 
is nowadays well founded. Most of the techniques used in the perturbative analysis 
to describe the parton evolution, adopt the tools developed by Dokshitzer et al. pf 



Amati et al. [27], Mueller [25], and numerous others (for an overwiew see ]29| , |30[]). 
The new component here is the extension to incorporate a space-time description on 
top of this formalism, which is commonly considered only in momentum space. 

It is evident that this paper attempts to join theoretical tools and concepts from rather 
different fields. Such a synthesis is necessarily a difficult task, and the the present initiative 
should be viewed as a first step in this direction. However, I believe that it is a promising 
approach towards a well founded and consistent description of the statistical properties of 
non-equilibrium parton systems. From the phenomenological perspective, it is an inevitable 
necessity to address this problem, since the experiments carried out at the HERA, RHIC 
and LHC accelerators will penetrate increasingly the physics of high-density QCD, where 
quark-gluon transport phenomena are of fundamental importance. 



2. FUNCTIONAL FORMALISM 

The aim is to describe the time evolution of a general non-equilibrium quantum system 
consisting of an ensemble of quarks and gluons in phase space, starting from some given 
inital state at time io- Since I am interested in the state of the system at finite times t > to, 
without a priori knowledge of the asymptotic final state at t = oo, the usual 5-matrix 
formalism of quantum field theory, based on in-out matrix elements, cannot be applied. For 
initial value problems as I want to describe here, the appropriate approach is provided by 
the functional integral formalism of the in-in generating functional for the Green functions 
of quarks and gluons, also referred to as closed-time-path (CTP) Green functions. The CPT 
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formalism is a powerful Green function formulation, originally introduced by Schwinger || 
and Keldysh j| for describing general non-equilibrium phenomena in field theory ||, ||, 
§]■ I 11 combination with the so-called non-local source theory and the loop expansion 
techniques developed by de Dominicis and Martin Q, and Cornwall, Jackiw, and Tomboulis 
[fLOf , one obtains generalized Dyson-Schwinger equations which incorporate the initial state 
correlations and provide a systematic treatment of the quantum correlations to any order in 
h. Furthermore it allows to describe phase-transition phenomena and dynamical symmetry 
breaking, issues that I will not address here, but which are of central interest when studying 
the confinement dynamics, as intended to be presented elsewhere. In this Section, I will 
first review the concept of the in-in generating functional and the effective action for the 
CTP Green functions, and then derive the dynamical equations of motion. For additional 
reading on these techniques I refer to the extensive review of Chou et al. @ and to the 
instructive work of Calzetta and Hu Jj]] . 

2.1 Preliminaries 

Starting point is the QCD Lagrangian given in terms of the gluon fields A£ and the 
quark fields ip, if) (which are vectors in flavor space, tp = (Vv ■ ■ •))> 



(i7/-^-M^K 3 



^ + u^n , (i) 



where F£ u = d^A v a — d v ' + g s f a b c A^A^. is the gluon field-strength tensor. The subscripts 
a, b, c label the color components of the gluon fields, and g s denotes the color charge related 
to a s = g 2 /(47r). The T a are the generators of the SU(3) color group, satisfying [T a ,T|,] = 
ifabcTc with the structure constants f a b c . The indices i,j label the color components of the 
quark fields and m = diag(m M ,mrf, . . .). I will in the following exploit the fact that high 
energies the quark current masses mf can be neglected, which corresponds to the chiral 
limit where they are exactly zero. 

In general the Lagrangian (|l]) must also include the Fadeev- Popov ghosts as independent 
field degrees of freedom. However, I will work exclusively in a class of ghost-free gauges, 
namely the so-called space-like axial gauges, which are defined by the gauge condition 

n^A^ix) = n-A a = , (2) 

where a constant four-vector in the x° — x s plane near the forward lightcone such that 
n 2 > 0. It may be parametrized, e.g., as = (a + 6, 0,0, a — 6), with the condition 
n 2 = Aab <s^ 1. The associated gauge-fixing term is denoted by a function ^ a {A^) which I 
take as 

U^) = - ^2 d x (n ■ A a )d\n ■ A a ) . (3) 
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Here a is the gauge parameter that specifies the type of axial gauge. In particular, I will 
henceforth set a = 1 which is known as the planar gauge. In contrast to covariant gauges 
where £ a {A) = -l/(2a)(d ■ A a ) 2 , the class of gauges (||) is well known to have a number 
of advantages [^6], First, the ghost fields decouple from the gluon field and drop out. 
Second, the so-called Gribov ambiguity is not present in this gauge. Third, the gluon 
propagator involves only the two physical transverse polarizations, which will simplify the 
analysis considerably. Furthermore, it allows for a rigorous resummation of the perturbative 
series at high energies in terms of the leading logarithmic contributions and consequently 
leads to a simple probabilistic description of the perturbative parton evolution within the 



(Modified) Leading Log approximation (MLLA) g|, IB @ in Q CD - 
The classical action corresponding to (|l|) is represented as 

JL^,^ = = / (0) [^] + i { ^\4>M + i {int) [A^^^} 

where 



(4) 



j(o)L4^ 



d 4 xd 4 y 



D 



d xd y j ip^x) [5^(x,y)J ipj(y) 

d 4 x | g s j^^A^x^x) + g s f abc [d fl A^ a (x)}A^ b (x)A" c (x) 
+ g*fabcfab>c'A^ b (x)A UyC (x)A£(x)A u c ,(x)} , (5) 



with the kernels of the free parts /^[A^ 1 ] and given by 

-l 



D 



ab 

(0) fiv 



(x,y) 



-i 



S ab S\x-y) 
Sij 5 4 (x -y)vy- d x 



(6) 



where the quark current masses are set to zero here and in the following. The operator 
is a generalized D'Alembertian containing the remnant of the gauge fixing term of ([!]), 
which for the gauge (p|) with a = 1 reads 



= 5' 



n^d v x + n"d£ 



n 



□ , 



(7) 



with D x = d x ■ d x , d£ = d/dx^. The inverses of (||) are the free gluon and quark Feynman 
propagators, i.e. the expectation values of the time-ordered products of the free fields 
-i(TA fJi (x)A v (y)) ( o) and -i(Tip(x)iJ)(y))^), 



D 



ab 



(0)A"A X >y) 



„-ih(x—y) r d/j, u (k) 

e ab „ 

k z + ie 

-ip.( X -y) , 1 

(2vr) 4 3 j-p + ie 



(2tt)4 
d 4 p 



n • k 



(8) 
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It is noteworhty that the form of the gluon propagator [] arises from the sum over the 
two transverse gluon polarizations, d fJiU (k) = J2s=i 2 £ m(^> s ) ' e iX^> s )> having the properties 



a 31] 



v U2 



dT{k) 



n u k 







0) 



meaning that only the two physical polarization states propagate, with e^k^ 1 = 0. For 
comparison, in the covariant Feynman gauge, d^ u = , dfi = 4, and k^d^ = k u ^ 0. 

In going over from the classical action (Q) to a quantum field formulation, the fields 
become Heisenberg operators. Let me introduce a compact notation for the different field 
degrees of freedom /: R 



4> f := (^,V,V') = (A^,^ u ,^ u ^ d ^ d ,...) 



f = g,u,u,d,d,.. 



(10) 



The state of the system may be characterized by the Heisenberg field operator &h(x), 
where &h = ^i/[<^/] an< ^ x = (^^O- Its time evolution is determined by the Hamiltonian 
H = #(°) + of the system (d t = d/dt), 



d t <S> H (x) = i [H, $ H (x)]_ . 



fill 



Defining t = to as the initial point for the time evolution of the system, the associated 
Heisenberg state vectors obey 



where 



Uj(t,t ) = Texp 



Uj(t,t ) \<t>(t ) ] 



( dt'd 3 x' J(x')<f> H (x') 
J to 



(12) 



(13) 



T denotes the usual time ordering operator, and the external source J is understood as 
a sum over sources for the various degrees of freedom. Note that the adjoint Uj(t, to) = 

exp i J* d^x' J(x')&h(x')\ involves an anti-temporal ordering TL In the absence of 
external sources, the state vectors are time independent: \<fi(t)) = \4>(to)). 

Upon switching from the Heisenberg picture to the interaction picture, the time evolu- 
tion of the corresponding interaction picture field <&i(x) is determined by the interaction 
Hamiltonian H mt alone, 



d t <S>i(x) 



(14) 



1 The apparent singularity of d M „(fc) at n ■ k ~ k + — must be dealt with the usual ie-prescription, or 
by taking the principal value. 

2 Since the quarks and antiquarks are treated as massless here, the different quark flavors are, with 
respect to the strong interaction, merely copies of each other. 



8 



where is related to the Heisenberg <3?# field by 

$/(z) = 5(t s t„) $h(x) S f (Mo) 
and evolves explicitly in time through 



S(t, to) = T exp 



dt'H^ int \t') 



to 



(15) 



(16) 



According to (|l2])-(|i"6|), at t = to the Heisenberg picture and the interaction picture 
coincide, 3>#(io,£) = $i(to,x). Hence, the interaction picture field $i(x) can be expanded 
at to i n terms of a Fock basis of free particle states, the in-basis, 

d 4 p 



= E / 7 ^-J(p°)(2n)5(p 2 ) Y,{ e ~***fM + e^afasj) (17) 

f=g,q,q J ^> 



(1) „(2) 



, n 



(oo) 



nn 



n 



(i), 



o: 



(18) 



where the a J- (a/) are the corresponding creation (destruction) operators for the particle 
types / = g,q,q with definite momentum pi and spin Sj, and the nV are the occupation 
numbers of the particle states, and 

a>f(Pi,Si)\0) = , = (n^ \ af(pi,Si)a^(pi,Si)\n^ ) . (19) 

Thus, a general multi-parton state \(j>) at time to is given by a superposition of such states, 

\</>(to)) = C(n (1) ,n (2) ,...,n (oo) ) | , n (2) , . . . , n (oo) ) , (20) 

with scalar coefficients C. Alternatively, the initial state of the system at to can be charac- 
terized by the statistical operator, or density matrix, 



(21) 



p(to) = I <Kh) > ( <Kh) I ( /5o )ij = < n w I p(to) I ) , 

which in the Heisenberg representation is time independent, but in the interaction picture 
evolves with time according to 

i \H (int \ p] , 



so that 



dtp 



p(t) = S\t,t )p(to)S(t ,t) , 



(22) 



(23) 



where S is defined by (16). For instance, a general density matrix that describes any form 
of a single-particle density distribution at to is 



p(h 



N exp 



E / * 



X 



d 3 p 



(2vr 



y^pF f {t ,x,p) a}{p,s)a f (p 



(24) 
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where denotes the hypersurface of the initial values and Ft is a c-number function related 
to the single-particle phase-space density of particle species / at x with four-momentum p, 
and N a normalization factor. 



2.2 The CTP generating functional 

After these preliminaries let me turn now to describe the time development of the multi- 
parton state from the initial state \<p m ) = \<p(to)), continously through finite intermediate 
times to < t < t^, to some final state \(p out ) = |</>(too)) in the remote future (see Fig. 1). 
In the usual S-matrix formalism of quantum field theory one calculates the in-vacuum to 
out-vacuum amplitude Z[J] = {0 tn \0 out ) j , and from this, physical quantities corresponding 
to in-out S"-matrix elements of certain operators, assuming that the Fock space of the 
asymptotic out-states is the same as for the in-states (Fig. la), as e.g. in scattering theory. 
In the present case, however, the system evolves forward through finite points of time, and 
so the asymptotic out-basis \(j) out ) is not known before the solution to the problem. There 
is an arrow of time, leading to an irreversible evolution. Moreover, in general \Q in ) 7^ |0 ), 
as for instance in the case of a phase transition or spontanous symmetry breaking where 
in- and out-vacua are of different nature. 

These problems can be overcome by using the CTP-formalism based on in-in rather 
than in-out matrix elements || |TT[| , but otherwise uses the familiar techniques of the path- 
integral method for quantizing the theory. The in-in generating functional is defined as the 
m-vacuum to in-vacuum amplitude Z[J,p] = TrJ2 v {0 m \ l p)j{ { p\p\0 tn )j, including possible 
initial state correlations represented by the density matrix p at to, and a sum over a complete 
set of states ip at too (Fig. lb). With reference to Appendix A, where the relevant concepts 
are reviewed and applied to the case of QCD, I merely state here the resulting path-integral 
representation for the in-in, or CTP generating functional. It is given by the following path 
integral representation in 2-point source approximation: 

Z P [J»,jJ,K^,k] = e ^ P U^3,K^,k] 

= J VA*V1>cP$ a exp i(l[A^ a ,^ a ] (25) 

+ tfAZ + j«$ a + Ti><* + \ a^a v p + 

where I introduced a shorthand notation for the integration over the space-time variables 
to be understood in the functional sense, 



J(j) = J <FxJ{x)(t){x) , <pK<p = J <Fx4)(x) K{x,y) (f){y) . (26) 

The CTP generating functional (25) differs from the usual generating functional of QCD 
in two essential aspects: 
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First, it contains both, local sources (J,j,j) and non-local 2-point sources (K,k). The 
former represent not only the usual external source contribution J ex t(x), but also the source 
term local source term K(x) for a possible dynamical background field present already at 
initial point t = to, that is, J = J ex t + K, and similarly for j,j. The non-local sources 
K(x,y),k(x,y), on the other hand, represent the 2-particle initial state correlations at 
t = to- Both these source contributions ^| stem from the general non-trivial density matrix 
p(to) that defines the initial ground state. In the usual field theory formulation both 
these source terms are absent. As a consequence, the connected generating functional 
Wp = —i In Zp in ( p5| ) gives both the non-local connected Green functions gluons and 
quarks, including initial state correlations (denoted by D^ u (x,y), respectively S(x,y)), as 
well as possible local mean fields which physically can arise either through non-vanishing 
external sources, or, in the case of gluons, may be generated dynamically by the system 
itself depending on the initial conditions (denoted in the following by A^x)). 

Second, the CTP functional Zp is defined on a closed-time path in the complex i-plane 
(indicated by the subscript P). This path P for the time integration is illustrated in Fig. 
2a: the path goes forward from to to too ° n the positive branch, and then back from too 
to to on the negative branch. Accordingly the generalized time-ordering Tp is defined such 
that any point on the negative branch is understood at a later instant than any point on 
the positive branch. This is not merely a mathematical trick to restore analogy with usual 
quantum field theory, but provides the means to compute expectation values for physical 
observables at finite time in contrast to the S-matrix formalism. The interpretation of this 
closed-time path is simple: although for physical observables the time values are on the 
positive branch, both positive and negative branches will come into play at intermediate 
steps in a self-consistent calculation, corresponding to a quantum mixing of positive and 
negative energy solutions. Therefore, in contrast to the usual path-integral formulation of 
quantum field theory, the non-local 2-point Green functions D fiv (x, y) and S(x, y) for gluons 
and quarks, respectively, come each in four different forms corresponding to the possible 

time orderings a(3 = ++, H — , — h, along the closed-time path P, as illustrated in Fig. 

2b. In as much as the propagators D flll {x, y) and S(x, y) can have values x and y on either 
the positive branch or the negative branch on the contour P, it is convenient to represent 
them 2x2 matrices G(x,y) = D^^S with components G a @ (a convention which holds for 
any 2-point function defined along the closed-time path P), 



3 In the 2-point sources approximation, the actually infinite series of non-local n-point sources that 
generate n-particle correlations, is truncated beyond n = 2 (c.f. Appendix A). 




(27) 
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where (...) = (cf) + , to\ . . . \<f>~, to) denotes the vacuum expectation value, if \(f>,to) = |0), 
or else the appropriate ensemble average. The generalized time-ordering operator Tp is 
defined as TpA(x)B(y) := 9(xo,yo)A(x)B(y) ± 9(yo, xo)B(y)A(x), where the +(— ) sign 
refers to boson (fermion) operators, and 6>(xo>2/o) = 1 (0) if xq > yo (xq < 0) on P. Hence, 
Tp coincides with the usual temporal ordering T on the positive branch (to —> too) of the 
closed time path in Fig. 2, but represents anti-temporal ordering T^ on the negative branch 
{too ~^ to)- The notation on the right hand side expresses that G F is the usual Feynman 
causal propagator, G F is the corresponding anti-causal propagator, and G > (G !< ) is the 
correlation function for xq > yo (xq < yo)- Explicitly, 



D ^ v {x,y) 



x,y) 



-i(TA^(x)A u (y) 
-i ( A^x) A v {y) ) 



D uu( x ^y) 



D ^(x,y) 



+i{A v {y)A ll (x) ) 
-t(Tt A IM (x)A u {y)) , (28) 



and 



S F (x,y) 
S < (x,y) 



-i ( Trp(x)ip(y) 
-i ( ip(x)tp(y) ) 



S > (x,y) 
S T (x,y) 



-i( tp(y)ip(x) ) 
-i(T^(x)^(y) ) . (29) 



The CTP generating functional Zp = exp(iWp), eq. (|2q), is the fundamental starting 
point for deriving the dynamical equations of motion for both gluon mean field A^, and 
the dressed gluon- and quark propagators, D^ v and S, using the matrix representations 
(f27|)-(p9"|). Formally, this Green function formalism on the closed-time path is completly 
analogous to usual quantum field theory, except that all propagators, self-energies, etc., 
are now 2x 2 matrices, as diagramatically represented in Fig. 3. Correspondingly, the 
Feynman rules remain the same, but each propagator line of a Feynman diagram can be 
either of the four components of the Green functions. 



2.2 The CTP effective action 

To proceed, it is convenient to work with the CTP effective action Tp, the tooparticle 
irreducible vertex functional, which determines the equations of motion for the physically 
relevant Green functions and the mean field, rather than with Zp or Wp of (25) which 
involve the sources J, K that do not have any immediate physical interpretation. The CTP 
effective action Tp is defined as the multiple Legendre transform of Wp |6|, [/J, which with 
respect to the 2-point source representation (p5t) is given by 



r P [^,r^,s] = w P [j»,j,j,K^,k] 



1 



A^K^ A" 



ipkip ) . 



(30) 



12 



Note that Tp reduces to the usual effective action for the one-particle irreducible vertex 
functions in the limit of vanishing mean field = and absence of initial state correlations 
K = k = 0. In the general case, one obtains 



r P [A fl ,D^,S] 



i Tr 



i 

-Tr 
2 



ln(£> ( -)£>) - D^D + 1 



(or 



HS7 ]S) - S7 ]S + 1 + T {2) [A^D^,S] . (31) 



The first term is of order h° and is given by the classical action (Q) (and eq. (|32| ) of 
Appendix A) with 



(32) 



The second and third terms are of order h 1 and correspond to the gluon and quark con- 
tributions in which the bare propagators and (||) are modified by the presence of 
a local gluon mean field A^ leading to 'mean field dressed' propagators and with 
an effective screening mass fl = ^[A^] (see Fig. 4). In analogy to (0): 



(D- ]r(x,y) - ^(x,y)5\x-y) 
S[ ](x,y) - jl q {x,y)8 i {x - y) . 



(33) 



Thus, the effect of the mean field is to shift the pole in the bare Green functions (||) by a 
dynamical mass function \x. 

The last term T^P in (|3l]) represents the sum of all two-particle irreducible graphs of or- 
der h 2 ,H 3 , . . . |u| ) with full propagators D^ u and S, dressed by both local mean field and 
non-local self-interactions (see Fig. 4). As will become clear, the real (dispersive) part of 

(2) 

T p contains the virtual loop corrections associated with the self- interactions of gluons and 
quarks, whereas the imaginary (dissipative) part contains the real emission, absorption, 

(2) 

and scattering processes. In other words, T P ' embodies all the interesting quantum dy- 
namics that is connected with renormalization group, entropy generation, dissipation, etc.. 
Explicitly writing out the color indices, it is given by (see Fig. 5a) 



Tr 



d zid z 2 \™u>% k v „ v , v {z 2 , zi;y) D% (x, z x ) D^„ (x, z 2 ) D^(y, x) 



+ 



d 4 zid 4 z 2 7 M 7$ Ejji u (z 2 , z^y) D%(x, z{) S^j^x, z 2 ) Sij(y, x) 



(34) 



Here A vu r u " and E. u are the qqg and qqg vertex functions, respectively, 



KX"{zi,z 2 ;y) 
E a iju {zi,z 2 ;y) 



\tSs\y-zi)s\y-z2)g[A»(y)] + o( g 2 s ) 

lv Tj 5\y - Zl )5\y - z 2 ) g[A»(y)] + 0(g 2 s ) , 



(35) 
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with A""/"// and j^T^ the corresponding bare vertices, and the function g[A^] describes 
the effect due to the presence of the gluon mean field as compared to free space, where 

g[o) = i- 

2.3 The self-consistent equations of motion 

The dynamical equations of motion for the gluon mean field and the gluon and quark 
Green functions in the absence of external sources are now as usual obtained from variation 
of the effective action Tp with respect to its variables, and setting the external sources to 
zero. Hence by functional differentiation of Tp ( |3T|) with respect to the gluon mean field 

one gets the Ginzburg- Landau equation J33| 



<5r P 



+ 2iTrJ 5[D <»> 



1-1 



l^(o)]- 1 



+ 



6T 



(2) 







(36) 



5A»(x) 8A^{x) I 5A»(x) 2 5A»(x) | 5A^{x) 

Similarly, the variation of Tp with respect to the dressed propagators D^ v and S gives the 
CTP version of the Dyson-Schwinger equations [34|, 



5T P 



SD^ u (y,x) 

. ST P 

— z 

SS(y,x) 



DfJ;{x,y) 



D 



(0) fiv 



(x,y) + U^ u (x,y) 



S 1 {x,y) - S^(x,y) + E(x,y) 



. 



(37) 
(38) 



Here n and £ are 2 x 2-matrices analogous to (27), representing the proper self-energy 
parts of gluons and quarks. They are obtained by functional differentiation of the quantum 
contribution T^P to the effective action (jsT]), 

,(2) 



2i dLp 



SD^x) 
From (04|), one gets (c.f. Fig. 5b), 



n£(*,y) 



5T\ 



SSij(y,x) 



T, i:j (x,y) . 



(39) 



-ids 



d'z l d i z 2 XZZ,A$, V X(»*,zi;v) (x,Zx)D^, (z 2 ,x) (40) 
+ J d 4 z 1 d 4 z 2 7 M T^, r^-, u (z 2 , z x ;y) Sij(x, z\) S V ji{x, z 2 ) 
Eij(x,y) = +ig 2 s Jd 4 z 1 d 4 z 2 j^T?,r b jjlu (z 2 ,z 1 ]y)D^(x,z 1 )S i > j/ (x,z 2 ). (41) 



The Dyson-Schwinger equations (|37]), ( pq ) can be brought into a more familiar form by 
employing the expressions for the free propagators <M) 



□ 



□ 



6a b g^S P (x,y) 
S a b9iJ,uSp(x,y) 



d"x'U^ aJ) ,(x,x')D^ b (x',y) 
d 4 x'D» a ^x,x')TLZ(x',yX42) 
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and 

j4_/ 



ij- d x - fiq(x,y)\ Sij(x,y) = 5ijd P (x,y) + J drx'E ik (x,x')S k j(x',y) 
Sij(x,y) -ry- d y - p, q {x,y) = 5ij5 P (x,y) + / (fx 1 S ik (x, x') T, kj (x', y) , (43) 



p 

where = d/dx^, nvf is defined by (0), the time integrations on the right hand sides 
are understood along the contour P, and the generalized <5p-function is defined on the 
closed-time path P (Fig. 2) as 

+<5 4 (x — y) if xo and yo from positive branch 
Sp(x,y) := ^ — S A (x — y) if xq and yo from negative branch • (44) 
otherwise 

Let me emphasize once more the essential difference to usual quantum field theory: the eqs. 



(42) and (43) are matrix equations and represent four equations, one for each of the four 
correlators (pij), respectively (|29|). In the limiting case where correlations among different 
partons vanish, one has G > = G < = 0, and because G F = G F t, one recovers the standard 
Dyson- Schwinger equations in terms of the Feynman propagators alone. The first equation 



in (42), respectively fl43|), describes the change of the propagators in the argument x, 
whereas the second equation describes the change in y of the adjoint propagators (adjoint 'f 
means hermitian conjugate with simultanous exchange of the arguments). A diagrammatic 
representation of these Dyson- Schwinger equations for the fully dressed Green functions 
Dfiu(x,y) and S(x,y) is shown in the previous Fig. 4. 

Let me summarize the considerations of this Section. The CTP generating funcional 
Zp involving a initial state correlations of the form (|2l]) , described by the density matrix p 
at t = to, yields infinite hierarchy of n-point Green functions, defined along the closed time 
path. As explained in Appendix A, the truncation of this hierarchy beyond n > 2 assumes 
that the dynamics may be described to sufficient accuracy by a possible local gluon mean- 
field and the non-local 2-point Green functions of gluons and quarks, and that higher order 
correlators are negligible. The resulting CTP effective action may then be represented by 
a systematic loop expansion corresponding to an expansion in powers of %. Considering 
the pure quantum regime with zero mean field, yields a coupled set of equations of motion 
for the gluon and quark propagators, which are 2x2 matrices containing the four possible 
time orderings of their arguments x and y. The solution of these dynamical equations then 
boils down to the evaluation of expectation values involving the propagators and vertex 
functions, e.g. by using perturbation theory M, |l0[|. 
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3. QUANTUM KINETIC THEORY 



Within the 2-point source approximation to the full theory in terms of 2-point Green 
functions, the resulting CTP Dyson- Schwinger equations (|4"2|), d43| ) ) contain the quantum 
dynamics in terms of the dressed gluon and quark propagators D^ v and S. Even with 
the neglect of higher-order correlators, the equations of motion are non-linear, non-local 
integrodifferential equations, generally not solvable in closed form. To make progress, one 
needs to supply reasonable physical input that allows to make realistic approximations for 
multi-parton systems of interest. 

First of all, I will confine myself for the remainder of the paper to the pure quantum 
dynamics of gluons and quarks, when a gluon mean field is absent. That is, I choose the 
homogenous initial condition A^(x) = at to, which in the absence of external sources 
implies that will remain zero at all times t > to, 

5T P 



5A»{x) 



0, A^ix) = (45) 



Consequently, in (plj), the classical part I[A^\ = (see also eq. (|32]) of Appendix A), and 
ft-g = ft-q = 0, so that the mean-field propagators reduce to the bare propagators, 

D^(x,y) = D^(x,y) , 5gjj(z,y) = S {0) (x,y) . (46) 

This step however is not an approximation, but merely serves as a simplification in order 
not to overburden the following analysis. The more general case including a dynamical 
gluon background field causes in principle no severe additional complexities, and will be 
addressed elsewhere. 

The essential approximation now is based on the 'two-scale' nature of high-energy QCD, 
as mentioned in the introduction. The dynamical evolution of a multi-parton system can 
- in a reference frame where the partons move with highly relativistic velocities - be char- 
acterized by two different time- (or length-) scales, separated by time dilation and Lorentz 
contraction effects: a quantum field theoretical scale Ar qua and a statistical-kinetic scale 
Ar^j n . This is illustrated in Fig. 6a. The quantum length scale Ar qua , measures the spa- 
tial range of quantum fluctuations, associated with the partons' self-interactions, and thus 
specifies the Compton wavelength A c = \i~ q of dressed partons. These gluon emission and 
absorption processes, embodied in the self-energy, dress up the bare propagators and allow 
to describe partons as quasi-particles with finite spatial extent, but with a dynamical sub- 
structure. This is nothing but the underlying philosophy of the usual parton description 
in QCD. The kinetic length scale Ar^, on the other hand, measures the range of binary 
interactions between these quasi-particles. These scattering processes may be described 
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on a semi-classical level, provided the local density of the quasi-particles is smaller than a 
critical density where the particles begin to overlap and the separation between quantum 
and classical regimes breaks down. Quantitatively one has to require that the mean free 
path X m f of particles is large compared to the radiative corrections to the Compton wave- 
length A c . The crucial point is, that with increasing energy scale the latter range becomes 
increasingly short-range, concentrated around the lightcone (see Fig 6b). Hence, in most 
physical situations at high energies, the quantum and the kinetic scales separate to very 
good approximation, and in the asymptotic limit exactly. It is important to stress that 
both quantum and kinetic scales define the microscopic regime of a semi-classical particle 
description. It is to be distinguished from the macroscopic domain of the dynamics of the 
bulk parton matter, characterized by comparably large space-time distances of the order 
n -1 / 3 , or n/(d r n), where n(r) is the density of quasi-particles. In this regime the system 
may be described by, e.g., hydro dynamical evolution, which is however beyond the scope 
of this paper. 

To exemplify this concept, consider the simple case of a parton in a Lorentz frame in 
which it moves with large momentum k + = E + k z (k +2 k 2 3> 1 GeV 2 ) nearly with the 
speed of light along the forward lightcone, x + = t + z. The quantum fluctuations around 
this parton's classical trajectory stem from its self-interaction with the gluon radiation 
field, corresponding to gluon emissions and reabsorptions, that smear out its energy over 
an interval AE ~ k 2 /k + . It may thus be pictured as an unstable particle with a typical 
life-time Ar p ~ 1/AE. On the one hand, in the direction parallel to the lightcone, the 
parton's intrinsic fluctuations decouple from the soft vacuum fluctuations with Ar„ ~ 
l/k + <C At p k + /k 2 fl35| . On the other hand, in transverse :r_|_-direction, the partonic 
fluctuations have a small spatial extent of Ar± ~ l/k± -C 1 GeV -1 . Therefore, on kinetic 
scales Akin > M > the parton appears as a dressed particle which can be considered quasi- 
classically as an extended object with a small transverse size Arj_ and a comparably long 
life-time Ar p — a quasi-particle. On quantum scales, however, the dressed parton has a 
substructure, determined by its surrounding cloud of gluons that it emits and reabsorbs 
due to its quantum nature. 

In this spirit I will classify the parton dynamics with respect to elementary and quasi- 
particle excitations, referring to them by the terms bare and dressed partons, respectively: 

(i) bare partons are to be understood as pointlike, massless quanta in the absence of 
radiative self-interactions, i.e. before renormalization. 

(ii) dressed partons, on the other hand, are dressed by the quantum self-interactions 
with their radiation field, which renormalize their masses and couplings. 

In the field-theoretical parton language, a dressed parton with its dynamically generated 
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renormalized mass can be described (in a frame where it moves close to the speed of light) 
as a bare quantum which is surrounded by a virtual cloud of other bare gluons and quark- 
antiquark pairs with which it emits and absorbs. Hence, kinetic space-time scales, a dressed 
parton can be visualized as a quasi-particle, i.e. an extended object with a dynamical 
substructure that is determined by the short-distance quantum fluctuations. 

3.1 Definition of quantum and kinetic space-time scales 

The realization of the two space-time scales, short-distance quantum and quasi-classical 
kinetic, allows to reformulate the quantum field theoretical problem as a relativistic many- 
body problem within kinetic theory. The key element is to establish the connection between 
the quantum-theoretical Green functions and the kinetic particle description in terms parton 
phase-space densities. In particular, the aim is to describe the evolution of a multi-parton 
ensemble, given at time to, with a certain spatial and a momentum distribution, by exploit- 
ing the 'two-scale' nature of high-energy QCD. As explained before, this requires a choice 
of Lorentz frame, in which the quanta move very fast and the typical momentum scale of 
their binary interactions and associated radiative processes is sufficiently large, such that 
the corresponding interaction times are small compared to the mean free time in between 
mutual collisions. For example: Imagine a high-energy reaction has produced an initial 
configuration of materialized partons (e.g. a hadronic or nuclear collision with y/s > 100 
GeV per hadron). If to denotes the earliest point of time in the lab frame, when the parton 
densities have evolved to satisfy ApAr 3> 1, where 1/Ap measures the scale of the partons' 
intrinsic quantum motion and Ar the space-time variation of the system of partons, then, 
for times t > to, an approximate incoherent treatment of quantum dynamics and kinetic 
evolution is justified, as has been shown by McLerran and Venugopalan (24]]. 

With this physical scenario in mind, now suppose, space-time is discretized into cells, 
with their size chosen intermediate between quantum and kinetic scales such that the sep- 
aration between the two scales is optimal ]?]]. Then the correlation between different cells 
will be negligible, and only when two space-time points corresponding to the arguments of 
the propagators or self-energies lie in the same cell, the 2-point correlation will contribute 
(as explained in Appendix B). Consequently, in a given cell, one can by construction neglect 
spatial inhomogenities of the local gluon and quark densities of the multi-parton system. 
Within each cell, one may therefore describe the short-distance quantum dynamics analo- 
gously as in vacuum or homogenous media, whereas inhomogenities of the spatial parton 
distribution and relaxation phenomena associated with binary collisions become apparent, 
as one moves from cell to cell. In continous space-time, corrections to this discretized 
picture can be taken into account by a systematic expansion in terms of gradients of the 
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spatial inhomogenities of the parton distributions (eq. ( |64[ ) below) . 

In order to quantify this concept, let me clearly specify quantum and kinetic domains 
with respect to the cellular space-time. It is important to realize that both quantum 
and kinetic scales are of dynamical, 'internal' nature, i.e. determined by the multi-parton 
evolution itself. The classification of the two scales only makes sense in the presence of 
self- and mutual interactions. However, the class of high-energy parton systems addressed 
here is characterized by two corresponding 'external' scales: first, the large energy scale of 
the reaction that produces the initial system of large momentum partons, and second, the 
initial local density of partons in phase-space that depends on the type of reaction. The 
first property implies a large characteristic momentum transfer q\ = (k\ — k^) 2 of scattering 
processes {k\k2 — ► k[k 2 ) and radiative processes (k\ — » k^k^)- The second property, on the 
other hand, is related to the mean free path A m / and mean free time r m j of partons in 
between subsequent scatterings. If the latter are large compared to the typical space-time 
extent l/q± of the the scattering- and radiation processes, then an incoherent treatment 
of the binary collisions among partons, and of the partons' propagation with associated 
quantum fluctuations, is applicable. This condition may be characterized by the invariant 
mass scale /i gq , defined such that 

q\ > a4 > A-2. (47) 

The parameter fi 2 q can be interpreted as defining the minimum virtuality of a dressed 
parton, or correspondingly, its maximum size, or Compton wavelength, A c = fJ-gg, such 
that the applicability condition of the parton description is ensured. Consequently, the size 
of each space-time cell must be chosen large enough that the spread of the dressed partons' 
intrinsic quantum motion is localized inside its four-dimensional volume, but smaller than 
the mean free path of dressed partons in between scatterings. Accordingly, I define the cell 
size Ar^ = Ar°A 3 r by 

ji" 4 < Ar^ ee ^\r) « K~ Q % D , (48) 

where Aqqd — 0.25 GeV is the QCD renormalization scale. For example, a cell size Ar <, 0.1 
fm allows to resolve particles with energy-momentum > 2 GeV. One can then characterize 
the kinetic space-time evolution of the system by a velocity profile of cells i, located around 
the points rf, with four-dimensional cell volume in its restframe 
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Each cell carries a total momentum 

i'"(>;) ■■ X /.;'/ (5U) 

and a total invariant virtuality (the incoherent sum of parton virtualities), 



Q 2 {n) := k 



2 

'j 

3=1 



(51) 

where the sums are over all dressed partons j inside the cell i, i.e. those that are during 
a time slice Ar° = contained within A 3 r = /i -3 around space-time point r{. The 
corresponding local four-flow velocity is (n) = P^/P°. This cellular space-time picture 
is illustrated in Fig. 7a. 

The validity of the above cell picture is controlled by the condition that the different 
scales are well separated: 

P +2 {n) » Q 2 (n) > M 2 ( n ) » A 2 QCD , (52) 



where P^ = (P + ,P-,P ± ), P ± = P ± P 3 , P± = JPf + P% with P+ (P~) the light- 
cone momentum (energy), P 2 = P + P~ — P"j_, and the normalization of a cell state \P) 
is (P\P') = 2P + (2tt) 3 5 3 (P - P'). On the basis of f)|) and in terms of these lightcone 
variables, the four-momentum of a parton j can be characterized by only two variables, 
namely, its lightcone momentum kj = XjP + with fraction Xj of the total cell momentum, 
and its off-shellness (invariant virtuality) k 2 = k^kj — k 2 j_ = M 2 (kj). Its lightcone energy 
is kj = (k 2 + kj ± )/kj~ ~ (k + 2 » k 2 > kj_), and one has therefore 

k 3 = (kf, k 2 ) = ( Xj P\ k 2 ) (2vr) (k 2 - M\ k]) ) = -±- 2 ^ dk) . 

(53) 

The requirement ( |52"| ) together with ( [48| ) hence translates to the parton level as kfj 2 S> 
fc| > /i 2 (rj), for all partons j within a given cell around r,- L . 

Q 2 (n) > k) > fi 2 gq > ^{n). (54) 

Since A c = fi~ x characterizes the maximum size of dressed partons, the ratio /i 4 (ri)//i^ g 
determines the minimum fraction of volume occupied by dressed partons in the cell. The 
quantum and kinetic space-time regions can now be defined as 

Ar qua = [P + -\ fi- 1 } , Ar km = \jT x , A Q ^ D ] . (55) 

For large P + and Q 2 , quantum and kinetic length scales are well separated and the 
parton phase-space densities Fj may be locally represented as a convolution of the kinetic, 
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statistical density of dressed partons A// of type / = g, q, with the quantum theoretical 
spectral density Vf of each dressed parton describing the intrinsic density of bare parton 
states as its quantum substructure, as depicted in Fig. 7b: 

F f (n,k) = F f ( n ,k + ,k 2 ) = N f (r i ,P + ,fi 2 gq ) ® V f (ri, k + ,k 2 ) , (56) 

where the convolution of the statistical density A// of dressed partons at the scale (i gq with 
the spectral density Vf, is defined as the average over the local space-time volume f2(rj) 
around r,- L of the densities, 

with n(n) ~ ^(n)- 4 , P + = P + (ri), k + = xP + , z = x/y (0 < z < 1), and 

A//(r, yP + ,fi 2 gq ) = ^^S(k 2 -^ 2 gq ) (58) 

/'=S,9 /'=S,5 

This ansatz describes the multi-parton system on the basis of treating each individial dressed 

parton as a composite particle of type / with a substructure of number of bare quanta 
f 

n J j of type / , weighted locally with the total number of dressed partons Nf in a space- 
time cell. The spectral density V characterizes the intrinsic structure of a dressed parton 
state, whereas the quasi-particle density describes the correlations and scatterings among 
those dressed partons. As will become clear later, the spectral densities Vf can indeed be 
identified with the QCD parton structure functions. The crucial quantities that control 
the cellular resolution in space-time of the partons' substructure are the characteristic cell 
size /i _1 (rj), and the minimum resolvable virtuality [x 2 q > fi 2 of dressed partons in the 
cell, or alternatively, the fractional space-time volume occupied, AQ/tt = /i 4 / ' {Ng q ^ gq ) A , 
that determines how dense a cell may populated without the partons overlapping. Hence, 
the validity of the kinetic approximation, based on the separation of quantum and kinetic 
scales, is controlled by the choice of these quantities, which need not be constant but rather 
may be taken as space-time dependent, i.e. variable from cell to cell chosen such that the 
resolution is optimal. A convenient choice would be, for instance, 

/*(»"*) - t^gq (60) 

which I will adopt in the following for lucidity, keeping in mind that /j, gq is not a free external 
parameter, but rather is to be understood as a dynamical, possibly space-time dependent 
quantity, which in principle should be determined self-consistently from screening effects. I 
will not address this latter issue here. 
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3.2 Wigner transformation and the kinetic equations of motion 



Let me proceed, referring to Appendix B for details, by introducing center-of-mass and 
relative coordinates of two space-time points x and y, 

r = - (x + y) , s = x - y , (61) 

in terms of which one can express any 2-point function W(x, y) = D^ u , S, H^u, as 

W(x,y) = w(r + ^,r-^ = W(r,s) , (62) 

and introduce its Wigner transform W(r, k) as [12] 

W(x,y) = J -0^e~ lk - s W(r,k) , W(r,k) = J d 4 se lh - s W(r,s) , (63) 

i.e. one Fourier-transforms with respect to the relative coordinate s being the canonical 
conjugate to the momentum k. In the cell picture of space-time, the coordinate r is the 
cell index that labels the kinetic space-time dependence 0(Ar&j„), whereas s measures the 
quantum space-time distance 0(Ar qua ), as illustrated in Fig. 8. In homogenous systems, 
such as the vacuum, translation invariance dictates that the dependence on r drops out 
entirely, and the Wigner transforms then coincide with the momentum-space Fourier trans- 
forms of the Green functions and self-energies. In general, spatial inhomogenities can be 
systematically accounted for by performing an expansion in terms of gradients d r = djdr^: 

W{r + s,s) ~ W(r,s) + s • d r W{r,s) + 0[(s ■ d r f] . (64) 

For quasi-homogenous, or moderately inhomogenous systems, such that s ■ d r W <C W, the 
correlations between different cells will be small so that the propagators and self-energies 
accordingly vary only slowly with r. One may then truncate the series fl64] ) after the second 
term, and convert the quantum field equations of motion (|42|) and ([43| ) into a set of kinetic 
equations by first performing the Wigner transformation ( |63| ) for all Green functions and 
self-energies, and then taking for ( |42|) and (|43| ) the sum and difference of the two adjoint 
equations in their transformed representation. 

This procedure (see Appendix B) yields two distinct equations for each of the Wigner 
transforms and S with rather different physical interpretations, which I will refer 
to as renormalization equation and transport equation, respectively. The renormalization 
equations are obtained as 

k 2 - \dl} D%(r,k) = -d^(k)8 ab i P + \ ( {II, D} + Y^ + 1 GZ { - ] (65) 
\ { 7 -p, S ij (r,p)} + = Sijlp - l - ( [ 7 -d r , + \ ( {£, + l -T\^ , 
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where d fiu (k) is given by eq. (||), <9 2 = d r -d r , and [A, 
The transport equations are found in the form 



AB-BA, {A,B} + = AB+BA. 



k.d r DZ(r,k) 
{7 • d r , Sij{r,p)} 



[n,£>]_ 
[7 • v , s] 



+ 

ab 



- Q h 

4 ^ ab 



1.1 



ib 

P,S]_ 



(66) 



Kl 



+ \& 
4 



In ( |65| ) and (66) the self-energies II and £ are explicitly given by (pop, (|4l|), and the 
operator functions Q and J^ 7 , on the right hand sides, which include the effects of spatial 
inhomogenities to first order in the gradient expansion ( |64|) , are ^ 



d A £,^S 



d*n%,d k x D 



->uu{+) 



d A s , <9?S 



{d£K,d r x D™} + -{d*K,d k x D™} 

(67) 

^ + ) = {a p A S,^} + -{a T A S,^5} + . (68) 

For completness, I note that the equations for quark Green functions can formally also be 
brought in a more familiar quadratic form, similar to the equations for the gluon Green 
functions, which exhibits the mass- and drift-term on the left hand side of the renormal- 
ization and transport equation, respectively: 



p ■ d r Sij(r,p) 



(7 • p + s ) Sij ip + - 



(7 • d r ) Sij lp 



{tf,s} 

E 2 , S 



1.1 



+ 4 



1.1 



where £ 



<%£, and 



B (±) 



5 ( (7 • p + s) (^H + ^(+)) ± (^H - ^(+)) (7 • p + £)) 

i f (7- S) (jrH + ± (jrH - ^(+)) (7. J 



(70) 



As will be seen in the following, the renormalization equations (|55|) express the nor- 
malization conditions imposed by unitarity and renormalization group due to the quantum 
self- interactions, and redefine the bare quanta in terms of renormalized quasi-particles. The 
transport equations ( |66| ) on the other hand describe the kinetic space-time evolution of the 
system of quasi-particles and their binary collisions. 

The kinetic approximation of trading the Green functions G(x,y) with their Wigner 

transforms G(r,p), means in the picture of cellular space-time, that inside a given cell 

carrying the space-time coordinate r = (t, r) as a label, G(r,p) equals the translation 

4 Note, dr = 9/<9fc M acts on a function f(r, k) as the derivative with respect to the space-time coordinate, 
whereas d£ = d/dk^ and dp = d/dp^ refer to the variation of four- momentum. 
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invariant Fourier transform G(p) of G(x — y), but outside of the cell it is zero. In another 
cell r', the Wigner function G(r',p') is determined by a different translation invariant G(p'). 
Hence, when looking at the short-distance quantum fluctuations within a given space-time 
cell around r = (r°,r), one may approximate the spatial distribution of partons as being 
homogenous and constant over the cell volume, and describe the short-range quantum 
dynamics in a tranlation-invariant manner. With the same accuracy of approximation, one 
can neglect in the quantum regime binary parton collisions, provided the mean- free-path 
Knf = {vgqFgq) -1 i n terms of the parton-parton cross-sections a gq and the local density 
F gq , is large compared to the spatial spread of the quantum fluctuations which is typically 
of the order of 1/yfp 2 . Hence, the essential requirement p 2 3> A m j can always be realized, 
if the particle energies are sufficiently large. 

On the basis of these considerations, I first study the quantum theoretical aspects em- 
bodied in the renormalization equations (|65|) to obtain the renormalized gluon and quark 
propagators, and from this determine the momentum dependence of the phase-space den- 
sities F gq , eq. (|56]) , associated with the variation of the parton structure functions. Subse- 
quently, I will investigate the transport theoretical aspects of the statistical kinetic dynam- 



ics, described by transport equation (36), which determines the space-time variation of the 
phase-space densities F gq in terms of renormalized, dressed partons. 

3.3 The 'physical representation' and strategy of solution 



Within the kinetic approximation, the goal is to obtain the best possible approximation 
to the complete propagators G = D^ V ,S, starting from the corresponding free Wigner 
transformed Green function. The 'free-field' solutions of the four types of correlators in 



(I 



9]), namely G F , G > , G < , G F , are in their most general form given by || 



1 



k 2 + ie 



2vri F( ) g (r,k)S(k 2 -p, 



(r,k) 
(r, k) 



D 



2ni (-d^(k)) [9(+k ) + F {0)g {r,k)\ 5(k 2 - fi 2 gq ) 

2iri (-V(fc)) [e(-k ) + F {0)g (r,k)] 5(k 2 - fi 2 gq ) 
1 



k 2 — ie 



2iriF {0)g (r,k)5(k 2 - p 2 gq ) 



where d fiu (k) is defined by (|8|), and 



+1 



+ 2vri F (0 ) q (r,p) 5{p 2 - p 2 gq ) 



'To* 

>(0)< 



7 • p + ie 
+ 2™ F (0) Jr,p)S(p 2 - p, 2 ) 



SUr,p) 



2m 



7 • p — ie 



1 - F (0) q (r,p)\ 5(p 2 - n 2 gq ) 

+ 27 " F (0)q(r,P)S(p 2 ~ (4q) 



(71) 



(72) 
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The scalar functions F( ) and F( ) g are the free-field analogues of ([56]) with the spectral 
densities Vf replaced by unity, F( ) f(r,k) = Mf (r,k) (g) 1, i.e. the phase-space densities 
of gluons and quarks that measure the number of non-interacting quanta in a phase-space 
element d^rd 4 p at a given time t = r°. Their presence is a direct consequence of the CTP 
formulation which incorporates initial state correlations due to a non-trivial density matrix 
p(to), eq. (21), corresponding to F(to,r,p) ^ 0, as opposed to the usual quantum field 



theory description, where p(to) = |0)(0| and F(to,r,p) vanishes. It is evident that in this 
latter case G > = G K = and G F = G F ^ at all times, so that the dynamics is describe by 
the Feynman propagators G F alone. 

More suitable for practical purposes, one may employ instead of the set G , G > , G < ,G F , 
an equivalent set of the retarded (advanced) propagators G R (G ) plus the correlation 
function G . The latter are directly connected with physical observable quantities, and are 
commonly referred to as physical representation ||. The functions G R , G A , G c are obtained 
via the relations 

G R = G F — G > = G < — G^ , G A = G F - G< = G> - G F , 

G c = G F + G F = G< + G> . (73) 

Because the fourth possible linear combination G F —G > —G < +G F is always identically zero, 
the three physical functions G R ,G A ,G C form a complete alternative set that eliminates 
the overdetermination of the set G F , G > , G < , G F . The 'free-field' forms of G R , G A and G c , 
corresponding to the ones of ( f7lj ) and (72) are 



D R W( r k ) - S R(A) (r v) - - 



(°)^ V ' ; k 2 ±iek (°) v,yj 7-p±iep 



D<o)»„(r,k) = —2m (—d fMV (k)) 1 + 2F (0) g (r, k) 8(k 2 - ^ gq ) (74) 



Sg(r,p) = -27ri( 7 -p) l-2F, 0) Jr,p) 5(p 2 - n 



where +(— ) in the denominators corresponds to the R {A). Generally speaking, the retarded 
and advanced functions characterize the quantum nature of parton states, whereas the 
correlation function describes the phase-space occupation of these states. 

The preceding relations (|7^) are generally valid for any 2-point function defined on the 
closed-time path, and hence apply to not only to the free-field case, but also to the full 
Green functions G = D„ V ,S, as well as to the self-energies £ = 11^, X. In matrix form, 
the correspondance between the representation (|27| ) in terms of G (£), and the physical 
representation denoted by G (£), is given by a unitary transformation U, being a 2x2 
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matrix with U mn = 1/ v2 for ran 
G = UGU- 1 = 




11,21,22 and Wi2 = -l/y/2: 
£ = U~ l £U = 




(75) 



where, in subtle contrast to (73), 

£ A = £ F + £> = -£< - £ F , 

£ c = £ F + £ J 



£ l 



£' 



£ F + £< 
- £> . 



£' 



£ l 



(76) 



The great advantage of this physical representation is that the dependence on the partons' 
phase-space densities F g and F q is essentially carried by the correlation functions G , 
whereas the dependence of the retarded and advanced functions, G R , G A , is weak. In the 
free-field case, this separation of correlations is exact, as is evident from (|74|), such that 
the retarded and advanced functions do not depend at all on F g , F q . In fact, even in the 
general case of interacting fields, this advantagous property becomes very suggestive when 
rewriting the renormalization and transport equations, (p5|) and (|66[), in generic form for 
the individual Green function components, 



{G-^G*-G*} = -2(V-k 2 ) (G R 



G J 



{5£,V} + + {r,<5G},(77) 



rt-i r c 

^(0) ' u 



-2i k-d r G 



c 



£ c , 5G 



+ 



5£, G 



c 



(78) 



where 



5G = ReG 
V = ImG 



1 



R , riA 

i ( G R - G A 



-[G n + G J 



5£ = Re£ 
T = Im£ 



R £ A 



(79) 



are the real and imaginary components of the retarded and advanced Green functions and 
self-energies, whereas 



G 



c 



£ 



c 



-(£< + £ ; 



(80) 



are the real correlation functions and corresponding self-energies. The physical significance 
of the (|77|) and (78) is the following: Eq. (77) determines the state of a dressed parton with 
respect to their virtual fluctuations and real emission (absorption) processes, corresponding 
to the real and imaginary parts of the retarded and advanced self-energies. Eq. (|7^), on the 
other hand characterizes the correlations mong different dressed parton states, and the self- 
energies appear here in two distinct ways. The first two terms on the right hand side account 
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for scatterings between quasi-particle states, i.e. dressed partons, whereas the last two terms 
incorporate the renormalization effects which result from the fact that the dressed partons 
between collisions do not behave as free particles, but change their dynamical structure 
due to virtual fluctuations, as well as real emission and absorption of quanta. For this 
reason £ R ( A ) are called radiative self-energies, and £ c is termed collisional self-energy. As 
shown by Kadanoff and Baym || , the imaginary parts of the retarded and advanced Green 
functions and self-energies (|79l) are just the spectral density V, giving the probability for 
finding an intermediate multi-particle state in the dressed parton, respectively the decay 
width T, describing the dissipation of the dressed parton. The general solution for V is 
given by 

= P - w+ WJW s APs£ + APr ' (81) 

where the second form in terms of the 'wavefunction'-renormalization (AVgs) due to virtual 
fluctuations, and the dissipative parts (AVr) due to real emission (absorption) processes, 
will prove convenient later ^| . The spectral density V satisfies the sum rule [||, || 

1 = Wl ^ T++V ^ k+ ^ = WI ^ k+ nr,k\k*) , (82) 

which is an implicit consequence of unitarity, and requires that the total lightcone mo- 
mentum of the spectral density of the internal bare partons must be equal to the dressed 
partons' momentum. For example, in the 'free-field' case, i.e. in the absence of interactions, 
one has AVse = <5(1 — k + /P + ) and AVr = with k 2 = fi 2 , so that V — ► V($) describes a 
single 'on-shell' parton state 

V {0) (r,k+,k 2 ) = d(k 2 - p? gq ) d( 1 - |L J , (83) 

on the 'mass-shell' k 2 = fj,g q , and carrying the total lightcone energy k + = P + ^ . This is 
nothing but the fact that the presence of a pole in the Green function means the presence 
of a particle, stable if it occurs for real k 2 , unstable if it occurs for complex k 2 , as in 
the Breit-Wigner formula (|8l|). The generalization of ( |33"| ) to the case of interactions, in 
which, as advocated before, a dressed parton may be visualized as a substructured particle 
with a fluctuating number of bare quanta intermediately present in its wavefunction, is 
straightforward. A dressed parton has now a "blurred" mass shell, because its internal 
excitations fluctuate due to virtual and real emission (absoprtion) processes of its bare 

5 This formula holds for both, space-like (k 2 < 0) and time-like (k 2 > 0) momenta. If is k 2 space-like 
then the imaginary part AVr vanishes, so that V is purely real. On the other hand, if k 2 is time-like then 
both AVse an d AVr contribute, and so V is complex. 

6 Note that for the choice (|6o|), the fraction z = x/y in the defining equation for V reduces to z = x. 
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daughter partons. The spectrum of these quantum excitations will have a finite extension 
around fi 2 q , described by the real part 5£ of the self-energy, with a width T, described by 
the imaginary part and being inversely proportional to the life-time of the particular parton 



state. Hence, one may write formally instead of (83), 

V(r,k) = 6^k 2 - M 2 (r,k + ,k 2 )^j , M 2 = SS - ^ , (84) 

where 58 and T are given in terms of the real and imaginary parts of the retarded and 
advanced self-energies (|79|). This representation serves to maintain the analogy with the 
free-field case, for which one has an immediate intuition. However, instead of a simple 
mass-shell condition, the argument of the <5-function now expresses a non-trivial functional 
dependence of the spectrum on k + ,k 2 and, in general, on space-time r. The solution of 
this implicit equation determines the spectral density V, which is subject of Sec. 3.4. 

Once the spectral density is known, the correlation function G c is given by the generic 
expression 

G c (r,k) = -2m [1 ± 2N(r, k)] ® V(r,k) = -2m [1 ± 2F(r,k)] 5 (k 2 - 7W 2 (r,fc)) , 

(85) 

where +(— ) is for gluons (quarks). It reduces to the free-field form, when V is replaced 
by V(o), eq. (H), so that = -2m [1 ± 2J\f] 5(k 2 - /i 2 gq ) = -2m [l ± 2F (0) ] 5{k 2 - 

fx 2 q ) becomes an 'on-shell' distribution, as in eqs. (0)- The theoretical basis for the 



previous, more physically motivated, ansatz ( J56| ) for the parton phase-space distributions 
F, becomes evident now: It is the logical generalization of the free-field forms fl74| ) to include 
renormalization effects and dissipation in terms of non-trivial spectral densities, or parton 
structure functions, which embody the underlying quantum dynamics. In this sense the 
Wigner functions F are the quantum kinetic extension of the classical particle phase-space 
distributions. 

Following this strategy, I will now proceed on the basis of the factorized ansatz ( jSlf ) 
for the gluon and quark densities Ff in terms of the quasi-particle distributions Mf with 
the spectral densities Vf, i.e. the presumption that the separation between quantum and 
kinetic scales allows a distinct treatment of the intrinsic quantum fluctuations of dressed 
partons and the kinetic correlations among them. In contrast to (74), the poles of the 
retarded and advanced Green functions are shifted by the real and imaginary parts of the 
self-energies £ R ( A \ and in the expression for the correlation functions, the <5-function is 
replaced by the spectral density V . Introducing the scalar functions for the fl for the gluon 
and £ for the quark self-energies through 

Kb = S ab {WW ~ 9^ k 2 ) II = 6ij p 2 E , (86) 
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instead of ( |74| ) one has now 

DK (A) (r,k) = ~ d ^ k) . S R ^(r,p) = . 1 ' V . (87) 

' fc2(l-n«W) P 2 (l_ 7 .pSii(A)) 

2tt» {-d^(k)) [l + 2Af g (r,k)] ® 7> fl (r,fc) 
2™ (-V(fc)) [1 + 2F 9 (r,fc)] ^A; 2 - ^ 2 (r,fc)) 
2vri ( 7 -p) [l-2jV a (r,p)] <g> P g (r,p) 

2iri ( 7 -p) [1 - 2F g (r,p)] 6 (p 2 - M 2 (r,p)) . (88) 



5 c (r,p) = - 



The two-step strategy that I will follow in the next Sections is then: 

1) In Sec. 3.4, the renormalization equations (^) will be solved for the retarded (advanced) 
Green functions 



ab 



D nU) ab (r,k) = (D- Q ] R(A) ) ab - (U m) ) 
(^(A))..(r,p) = (*S ( o)V)L - ( S *w)« ' 



(89) 



which determine the spectral densities V g and V q in terms of the radiative self-energies 
Il^andX^). 

2) In Sec. 3.5 the transport equations will be solved for the correlation functions 



Sdj{r,p) 



u Raa' 



S f 



K.io)::-(no)^ 



U A b'b 



S (0)Cj i , jl ~ ( s c)i'j' 



S A 



(90) 



which determine the parton phase-space distributions F„ and F„ by the collisional self- 



energies 11^ and £ , in conjunction with the spectral densities V g and V q . 

3.4 Quantum dynamics and renormalization equations 

As advocated above, when addressing renormalization effects and dissipative quantum 
dynamics, it is appropriate to focus on the retarded and advanced propagators and the 
imaginary parts of the self-energies, which embody the short-distance propagation of quan- 
tum fluctuations. Furthermore, on quantum scales, one can neglect the r-dependence, and 
thus ignore in this regime the functions Q and J- in eqs. ( |65| ) and (|66|). Then, by perform- 
ing the transformation of the Wigner transformed Green functions D^ u , S and self-energies 



!!„!,, E to the physical representation via (|75[), one obtains simplified equations for the 
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retarded (advanced) functions (c.f. Appendix B). The renormalization equations reduce to 
the following form: 



k' 

4 



9r) D^ R{A) (r,k) = -d^{k)5 ab + X - (^U R ^D R ^ + D R W U R ^ 
\ {^P,S« {A \r, P )} + = Sa + \ (V^S*^ + S R Wx R Wy (91) 

To solve the equations (pl|), it is suggestive in view of fl74|), to parametrize the renormalized, 
'dressed' propagators D R , S R on account of their Lorentz structure as |26|, 

ET R(A) (r k) - 6 h -d» v (k) = A (r k 2 kA 5 fc Z^eM + 

D * (r ' fe) " dab kHl-fl R ^)(r,k)) ~ Afl(r ' fc > Kk)dab **±ieho + - - 

S$ A \r,p) = 5, ^ 



p 2 (l -7-pS R (^)(r,p)) 
= A q (r,p 2 ,K p ) 8ij P — + A q (r,p 2 ,K p ) 5ij?—^ + ... , (92) 

where d /ll/ (k) = g^ u — (n^k^ + n v k^)/{n ■ k) as before, the scalar self-energy functions Ii R , 
t R are defined by (pq), and /c implicitly accounts for the dependence of A 9 and A g on the 
coordinate r, which is conjugate to k + . The function k is of the order of the large lightcone 
momentum k + squared (c.f. Appendix C) 

nl = ^ ' ^ ~ k +2 , n 2 « 1 . (93) 

The renormalization functions A g (A g ) account for the modifications of the 'bare' propa- 
gators (74) due to the self-interactions embodied in II (£). The third function A„ turns 
out to be proportional to A q (c.f. Appendix C). They are normalized in accord with the 
condition (M), such that 



A 9 (r,k 2 , Kp ) = A q (r,p 2 , Kp ) =1, (94) 

meaning that a gluon or quark is considered as a maximally dressed particle (in the sense 

of the applicability of (54)), corresponding to the invariant scale /i gq , which may be called 

the dressed partons mass-shell. 

It is well known p(| that other contributions to the propagators in (92) , indicated by the 

dots, are strongly suppressed f\ in lightcone-dominated processes Q 2 —* oo by inverse powers 

of Q 2 . In fact, this feature is the very foundation of the QCD parton description within the 

7 Note the benefit of the employed gauge n^A^ = for the gluon field, eq. by suitable choice of 

the vector such that n 2 < 1, one can concentrate the short-distance quantum fluctuations to arbitrary 
proximity around the lightcone, n 2 — > 0, i.e. Kk — > oo, corresponding to the asymptotic linit k + — > oo. In 
this regime the leading log singularities of the propagators give the dominant contributions and the dotted 
terms in and (J9q) can be neglected because they do not generate leading logs. 
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(modified) leading log approximation (MLLA) p6| , p^ ] , where the renormalization problem 
reduces to a multiplication of the bare propagators and vertices by scalar functions. That 
is, as depicted in Fig. 9, with respect to the Wigner transforms of the self-energies (40) 



and (41), 



n 



livR{A) 
ab 
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d A k! 
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(95) 



(2vr) 4 i 



7^n S%f> (r, —k' + k) Z b nll , u (-k' + k, k, k') S R(A) (r, k') 



d 4 k> 



lc TZ v S R{A) (r, p - k') Ej,j T (p-k',P, k') D a J d H ^> (r, k') 



or R(A) . 



(2vr) 4 i 

one can represent the qqg- and ggg vertex functions T, respectively A, eq. (^), as multi- 
plicative renormalization functions, 



rf 7 a (r; Pi,p 2 ,k 3 ) = ^ T t a j V qqg (r; p 1 ,p 2 ,k 3 ,K 1 ,K 2 ,K 3 ) + ... 



tf%?\(r; h,k2,h) = f abc g^\ K V g99 (r; k{ , k 2 , h%, «i, «2, «*) + 



eabc 



2 ,2 ,2 



(96) 



where k v is understood as the momentum associated with the intermediate of the three 
gluon virtualities, i.e. k v = fav if k\ < k\ < fcf, etc.. Analogous to (|9^), the normalization 
conditions are 



Vqqg( r ', Pit K i 



1 



(97) 



Employing these definitions in the renormalization equations (|9l| ) and the expressions 
for self-energies given by ( |40| ) and (0^), and differentiating the inverse of the propagators 
D R y A ^ and S R<yA ^ with respect to the gluon and quark virtuality, respectively, 



J^(-dr(k) [D R ( A \r,k)]~l 
^(7-p[^)(r,p)]- 1 



5 



n^ A )(r,fc) 



2 t R{A \r,p) 



(98) 



one obtains the following determining equations for the momentum dependence of the scalar 
self-energy functions (|86|), fl R ( A ) and Y, R ( A ^ to order g 2 s in terms of the renormalization 
functions A g and A q , 

,2 , d Ay ,ro+l/(2 M ) 

; — \T~ / drTtk'r) 

{27r)H Jr -l/(2 M (r)) 

, mV .-; A; 2 , A;' 2 , fc" 2 , « fc , « fc ,, Kfc „) A fl (r, fc")) x k" , n) 

d 4 p ' /T +l/(2|») 



dA g (r,k') d 



dk 
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drTip'r) 



(99) 
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dA q (r,p') d 
dp' 2 dp 



TyTf [ V gqq( r > k2 iP 2 iP" 2 ^k,Kp',K p >') A q (r,p")^j xU% q (p',p",n) 

f) r r/ 4 ™' /•r + l/(2a) 

f*t«W(r,p) = -9 2 s C F (2*i) drT(p'r) (100) 

Op z J (27T) 4 Z Jr Q -l/{2fi) 

( V q 2 qg (r; p 2 ,p' 2 , k" 2 , K p , n pl ,K k n) A g (r, k")) x U™(p', k" , n) , 



dA q (r,p') d 
x " 12 ay/ 2 



dp' 2 dk 

plus terms 0(n 2 ) which can be neglected for n 2 -C 1. The constants CA$ab = facdfbcd = 
N c 5 ab , T R 5 ab = Tr(T a ■ T b ) = \5 ah , C F 5ij = (T a ■ T a )ij = %^<%, arise from summing 
over color indices, and Nf is the number of quark flavors. The time integral / dr on 
the right hand sides extends over the finite time slice fi^ 1 = fj, (r) of the space-time 
cell around r, weighted by the function T which incorporates the relation between 
momentum and space-time as constrained by the uncertainty principle: it limits the range 
of virtualities k' 2 such that within the finite time interval Aro = only those fluctuations 
k — > k'k" are resolvable that are sufficiently short-living, with proper life-time to — 1 jk' and 
7To — k + /k'2 < Aro. Finally, the functions Uj (k',k",n) represent the squared matrix 
elements for the virtual decay processes k — ► k'k" . 

In the cellular space-time picture, the momenta of partons in a given cell around r are 
per design limited by the condition (|54|), such that 



P +2 {r) > k +2 > k 2 > i? gq > fi 2 {r) . (101) 

As explained in Appendix C, by employing this condition and introducing the fractional 
lightcone momenta of the daughter partons in the process k — ► k'k" , 

k! + k" + 

zy = z = — , z k n = 1 - z = -p-- , (102) 

d 4 k' = - dk' 2 dk" 2 dze(k 2 -^-- . (103) 



the integrals (|9^) and ( |10C| ) are readily evaluated to leading log accuracy. Solving for the 
renormalization functions A g , A q , and taking K k = k +2 ,K p = p +2 from (|93[), the result is 
(c.f. Appendix C): 

A g (r,k 2 ,k+ 2 ) = exp|-£ + ^lj\ zA (r,z,k' 2 ) Q 7 fM + ->«(*, e))j 
A q (r,p 2 ,p +2 ) = exp(- /f %i f ^A(r,z,p' 2 ) 7 f(^e)) , (104) 



where 



^(r,A; 2 ,z) * " s((1 2 ; Z)fc2) - M (r)) (105) 



32 



is the effective local coupling strength 8 averaged over a cell centered around r of space- 
time extent O(r) ~ A 4 r, and a s {q 2 ) = b In (^Mqcz))] -1 ' 6 = 1UV C - 2N f . The function 
A(r,z,k 2 ) reflects the fact that in cellular space-time the relevant quantum fluctuations 
are restricted by the uncertainty principle, as embodied in the ^-function term: k 2 < [i 2 = 
l/(Ar) 2 , A; < 1/Ar°, with Ar° = k+/k 2 . 

The functions 7(2, e) in (104) involve, at leading-log level, the standard DGLAP kernels 
|p6| , carrying an explicit /c + -dependence arising from the dependence on k ~ k + 2 , 

z 




ifM = 

7 q g q (z,e) = 

7 f(z,e) = (', ( , •,',:,■■, ) !1(, <>> 



1 - z + e(/cjfc) 
where in the denominators the function e appears, 



£(K ' ) = a ' (107) 



which arises here as a consequence of the d 2 term in the renormalization equations (91), 
after Fourier transfoming with respect to r~ = r° — r 3 being the conjugate variable of k + . 
It can be interpreted as manifestation of the indeterminancy principle, which determines 
space-time uncertainty of the order of the cell size Ar that is associated with the off-shellness 
of the partons. The presence of e effectively cuts off small-angle gluon emission when the 
emitted gluon is soft, i.e. when z g = 1 — z — > 0, by modifying the free gluon propagator 
oc 1/zg to the form l/(z g + e) when k/k + = O(l), that is, in branching processes with large 
space-time uncertainty. This ensures that the two daughter partons can be resolved as 
individual quanta only if they are separated sufficiently by Ar cx 1/k in position space, in 
accord with the uncertainty principle. Note that e can be neglected in the terms oc l/(z g +e) 
in ( |106D for energetic gluon emission (z g — ► 1), but is essential in the soft regime (z g — > 0). 
The effect of e has been shown [26, to result in a natural regularization of the infra-red- 
divergent behaviour of the branching kernels ( |106| ), due to destructive gluon interference 
which becomes complete in the limit z g — > 0. 

As summarized in Appendix C, the renormalization functions Ay are intimatly related 
to the spectral densities V g and V q defined by ( ]59| ) and (pi]), 

f'=9,Q 



8 In (105) one could imagine, instead of the ^-function 6(k 2 — (J 2 ), a smeared-out probability distribution, 
e.g. oc exp(— k 2 /fi 2 ), by choosing a more refined form for the function T under the time integral in eqs. ( p9| ) 
and (100). The specific choice is ambigous at this level of calculation (see Ref. |36| for details). 
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v * = E *f = * Tr 



(108) 



where the Vi have the following solution, 



V f f '(r;x,k 2 ) = S f f ' 5(1 - x) 5 (k 2 - fi? gq ) A f (r-p? gq ,k +2 ) 



+ Af(r,n 2 gq ,k +2 ) E/^ W /> \A(r,k 2 ,z) 7 ff'(z,e)vf, (p*' 2 5*' +2 



f" ^99 



x AJ 1 L^-,k' + 2 



(109) 



Comparing this expression with eq. (|8l|), V = 5(k 2 — {j, 2 )AV$g + AVr, one sees now that the 
Sudakov formfactor A f together with the real emsission and absorption probabilities Wf = 
J2f" I dink dzA^fj Vet, combine to play the role of the 'wavefunction renormalization' 
part AVse an d the dissipative part AVr- Equation ( |109| ) has a simple physical significance: 
The first term on the right hand side represents the probability to find a dressed parton 
of type / in the cell around space-time point r as a 'classical' particle, i.e., without any 
other gluons or quarks present in its wavefunction, or spectral density. In accord with the 



normalization (94), this means that it is propagating 'on-shell' with k = [i > jj, (r) , where 



^(r) _1 is the resolution size of the cell as explained in Sec. 3.1. Its fraction of the cell's 

lightcone energy is x = k + /P + (r) = 1. The probability for finding such a rare fluctuation 

is suppressed by the functions Aj(r,k 2 ,k +2 ), which becomes stronger with increasing gap 

between fi 2 q and k 2 . The second term on the right hand side corresponds then to the adjoint 

probability, that the parton is actually a dressed parton with a substructure, described by 

the balance between real and virtual emission and absorption processes, while localized 

within the cell around r. It is obvious, that the spectral densities of dressed partons, 

introduced Sec. 3.3, are identical to the usual parton structure functions Vf, i.e. the 

probability densities for finding a dressed parton / in an intermediate state containing 

number of bare partons as virtual and real fluctuations. 

From ( |109|) with (|104| ), and using the representation (|56| ) of the parton densities Ff in 

f 

terms of the parton structure functions Vj , follows then the final form of the renormaliza- 
tion equations, 



*wM^) = /„' 



dz A(r; k 2 , z) 



-F n 



r; — , zk 2 

z 



Ik 



r; x, k 2 ^j 



+ 



2N f F q (r;x,k 2 )T^(z,e) - N f F g (r;x,k 2 ) Tf(z,e) 



dp 



x,p 



dz A(r;p 2 , z) 



1 r, ( x , 
-F q \ r;-,zp' 

z \ z 



F, 



(r;x,p 2 ^ 



Tf(z,e) 
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+ F g (r;x,p 2 )rf(z,e) } , (110) 



which are the space-time generalization of the DGLAP evolution equations [26] that govern 
the momentum dependence of the parton densities. The effective branching kernels Ij J 
[j37f are related to the by, 

where the upper (lower) sign stands in the term in brackets is for gluons (quarks). It yields 
a suppression, when the phase-space density F g or F q becomes large, so that the emission 
processes / — > /'/" are competed significantly by absorption processes /'/" — ► f. In the 
limit Ff/ 3> 1, the detailed balance is established, in accord with the Bose-Einstein and 
Fermi-Dirac statistics of the gluon, respectively quark densities. For instance, in thermal 
equilibrium, F^T^ = {e~ E l T =F l) _1 i so that T = 7(1 — e E ^ T ), which tends to zero as the 
temperature T becomes large. 



Eqs. ( 110 ) are the main result of this Section. They emerge as a direct consequence of 
the renormalization equations (|6^) in the short-distance regime of virtual and dissipative 
quantum fluctuations, and ensures unitarity conservation locally in each space-time cell. 
They embody the Heisenberg uncertainty principle, expressing the fact that it is impossible 
to localize soft partons in a given cell if their wavelength exceeds the cell size, which sets 
the resolution scale. Last not least, they account for the balance between real emission and 
absorption processes that, in tends to increase (decrease) the effective real emission rate of 
gluons (quarks). 

3.5 Kinetic dynamics and transport equations 

With the dynamical structure of dressed partons quantitatively controlled by the above 



renormalization equations ( 110 ), one is now in the position to address the kinetic space- 
time evolution of the multi-particle system in terms of statistical binary scatterings among 
these dressed partons. As explained in Sec. 3.1, in order to obtain quasi-classical trans- 
port equations for the phase-space distribution functions F = N <S> V, two key conditions 
have to be met: First, as before, the maximal space-time extension of relevant quantum 
fluctuations, A c = \i ~ , , is supposed to be smaller than the mean free path A m / between 
scatterings. Second, the typical four-momentum transfer q± = \/\q 2 \ in the scattering of 
any two partons, is required to be larger than the inverse Compton wavelength A" 1 = \x gq . 



That is (c.f. eq. ©) 



Am/ > Vgq , Vl > V 2 gq = K 2 ■ ( 112 ) 
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The first condition ensures that the quantum evolution, taken care of by the renormalization 
equations, can be factorized from the scattering processes. The second condition guarantees 
that the scattering is sufficiently hard, i.e., is of short range compared to the space-time 
extent of the scattering partons' intrinsic quantum motion, so that over the duration of the 
scattering, the dressed partons can be treated as 'frozen' assemblies of bare particles that 
represent their instantanous quantum state (the usual sudden approximation). These two 
conditions are equivalent to the factorization assumption of the well established 'QCD hard 
scattering picture' []38| for, e.g., high energy hadron-hadron collisions, where the colliding 
hadrons are described as conglomerates of bare partons in terms of their structure functions. 
The relation of this hard scattering picture to the present approach is its adoption to 
multiple, internetted scattering processes in a system of stochastically colliding dressed 
partons, each of them represented by their own structure function, or spectral density. 

The two requirements ( |112| ) are the crucial points, which allow in the following to cast 
the kinetic evolution into simple, probabilistic Boltzmann-type equations, which however 
have to be solved self-consistently in conjunction with the renormalization equations ( |110|) . 
Here is the key difference from other formulations @, [3^, 40 1 to include quantum effects 
in a quasiclassical treatment of transport phenomena, in which one has only one type of 
equation, a generalized Boltzmann equation, that contains local, classical part and a non- 
local quantum contribution, containing the space-time history of memory effects. In the 
present approach, this is translated to stochastically occurring ('local'), hard parton-parton 
scatterings, linked with the causal quantum evolution in between scatterings ('non-local'), 
accounting for renormalization and dissipation due to the previously occured scatterings. 
The advantage is here, that while quantum effects are included in the multi-particle evolu- 
tion, still a local (in space and time) picture can be maintained, where memory effects are 
embodied effectively in the dressed partons' structure function evolution. 

To proceed, recall from Sec. 3.3, that the correlation functions D^ U ,S C are the quan- 
tities which are determined by the transport equations of the form ( f78| ) (explicitly given 
in Appendix B, eqs. ( |167|) and (|16S| ). On account of the presumed conditions (|112|) , over 
kinetic space-time scales X m f > A c , the quantum motion decouples, so that the correlation 
functions are determined by the collsional self-energies n^E^, in conjunction with the 
real parts of the retarded and advanced functions (c.f.(|78|)). This means that the collisional 
self-energies are to be evaluated with the renormalized propagators ( |9^ ) and vertices 
which were obtained from the retarded and advanced self-energies before. 

Noting that from (|73|), one has 



D° v = D> v + D< S c = S> + S<, (113) 
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the transport equations (|66|) now read (c.f. Appendix B), 



k-d r D% € (r,k) = - - (l^D A + U R D* - D^U A - D R U^\ (114) 



2 



(IV 



ab 



i { 1 -d r ,Sf j {r,p)} + = - ( [l-P,S^]_y + ^S A + Z R S* - S*Z A - S R ^ 



which can be rewritten as 



(IV 



ab 



k-d r D%*(r,k) = - l - nn>(r,k),D<(r,k)} + - { n<(r, k) , D>(r, fe)| 

p-drSf^p) = \ ^|s > (r,p),S' < (r,p)} + - |s<(r^),5 > (r,p)}J .(115) 

These equations correspond to what is usually termed the quasi-particle approximation. 
In the cellular space-time picture, the characteristics of the statistical kinetic evolution of 
the system are, per design, insensitive to the localized fluctuations associated with short- 
distance quantum dynamics inside a space-time cell. To stress it more precisely, the space- 
time variation can be considered homogenous over the range of the Compton wavelength 
A c = < /i -1 (r), so that G c (r, k)\ 3> \X^dfG G (r, p)\ and the derivatives with respect 
to r on the left hand side of the original eqs. (|65|) and ( |66|) may be omitted. In the 
present context, it emerges as the logical consequence, that the partons can be described on 
kinetic space-time scales as quasi-particles, with the underlying quantum motion effectively 
accounted for in the renormalized propagators and vertices. 

The self energies 11^, and are obtained from the general expressions ( fl"0| ) and (|4l"D, 
respectively. The lowest order non-vanishing contributions are the two-loop diagrams shown 
in Fig. 10, which are proportional to h and 0(g 4 ). In terms of the renormalized correlators 
and one finds: 



■ > 



g A s f d 4 k' d 4 q 



faa'c^(i(i'<r( k,k k ,k ) DZZ < (r,k ) 



2 J (2vr) 4 i (2vr) 4 i Jaac ' 

x fc'feK'p\{-k', -q, q + k') D p f p fl ^(r, -q) / e '/''d'A A y r /(-g - k', q, k') D T ^(r, q + k') 
x fdb'bKM-k', -k + k', k) D££*(r, k - k') 

+ l / (2^ <l',k- Q - q', q) DSHr, g) D5£*{r, </) (116) 

x vg&{-q, -k + q + q', -q', k) D% k {r,k-q- q>) 
+9t N f J ■ { ^)H /-tfcW-fc, k - k' , k') ug*(r, k') 1( ,,Ti S* n ,(r, -q) 

x 7r ,T n d ,V ^(r, q + k') D T J T ^(r, k') f db r b Kv'v(-k', -k + k' , k) D^\r, k - k') 
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d 4 k' d 4 q 
(2vr) 4 i (2ir)H 



^T? i S*Jr,k-k>)i v lf n S}, i (r,k>) 



> . 



x 7r T4, S< n , (r, k' - q) la T^, S< (r, k') D% * (r, q) 



9 



4 ' d 4 k' d 4 q 



la T^D^'\r,k') 



2 J (2vr) 4 i (2vr) 4 i 

x fc'feK'p\{-k', -q, q + k') D p f p f ,^(r, -q) f e 'f"d'^y P 'r'{-q - k', q, k') 
xD^(r,q + k') lT T^Sl,(r,p-k') 
d 4 k' d 4 q 



+9. 



(2vr) 4 i (2vr) 4 i 

x Sl n ,{r, -q) la ,Ti 4(r, q + k f ) lr ,T* v D^ d \r, k') 
x lr T^S^,(r,p-k') 

d 4 k' d 4 q 
(2tt)H (2vr) 4 i 



la nS? v (r,p-k')Diy<(r,q) 



x lxT§ m Sl m , (r, p-k'-q) 7t Z* D% * (r, k') 



where \ wv \p\1P21Vz) and v^^J (pi,P2,P3,P4) are the usual 3-gluon and 4-gluon vertices (c.f. 
Appendix D, eq. (|223j )). The correlation functions D^ u , are related to the phase-space 
densities F g ,F q by fl88|), i.e., 



D^(r,k) 



-2m (-d^k)) 



0(±k°) + F g (r,k) 



5(k 2 - M 2 g (r,k)^j 



-2m (7 • p) 



6(±p°) - F q (r,p)} 8{p 2 - M 2 g (r,p)) , (118) 



where the signs +(— ) refer to > (<) []. As repeatedly stressed, the densities F are the distri- 
butions of dressed partons, with their substructure represented in terms of corresponding 



assemblies of bare partons, that satisfy the condition (54). Therefore, the functions F 

can also be interpreted to measure the number of bare partons with dynamical invariant 

masses k 2 > fi 2 q . Consequently, a binary collision of two dressed partons can be described 

in terms of the above 'hard scattering picture' as a scattering of two bare partons, one out of 

each assembly, picked statistically from the instantanous quantum state of the two dressed 

partons, as given by their structure functions, or spectral densities. The four-momentum 

transfer q 2 ^_ sets hereby the probing scale, so that k 2 ~ q 2 ^_. Therefore the energy spectra of 

9 It must be mentioned that eq. ( [L1S| ) assumes a spin-symmetric form for the quark-antiquark spinor 
products, which means a neglect of spin-polarization effects. As shown by Elze et al. in general the 



quark phase-space distribution does require at least an 



x 



-matrix representation. 



(117) 
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partons emerging from these scatterings are - using eq. - determined by 



n, x ( M 2 Jr,k + ,k 2 )' 

±k°(r,k) = u (0) 1 + §y 



2 uj 



(0) 



# = ±P (r,p) = E {0) 1 + 



2 B 



(0) 



(119) 



where A: = l/2(k + + k~) ~ fe+/2, p° ~ p + /2, and w (0) = ±^k 2 + /i 2 gq , E {0) 
±^Jp 2 + /x^g. Accordingly, one can now write 



±Vfc 2 + «i 



I sjp 2 ■ T " 



2 \ 
H-gq) ) 



(120) 



and 



-Fs( r > A;) 5(Ar — q±) = — 
F q (r,p)6(p 2 - q\) = — 



Fg(r,k)5{k° -u) + F g (r,-k)5{k° +uj) 
F q (r,p) 5(p° - E) + FJr, -p) 5(p° + E) 



(121) 



which exhibits explicitly the particle-antiparticle character of the phase-space densities. In 
particular, F q (r,p) denotes the quark distribution and F q [r, —p) the antiquark distribution. 

Using the representations (|11S| ) for D^ u , in the self-energies (117) and ( |116|) , and 
substituting into the equations ( p.l5|) , gives the final form of the transport equations in the 
kinetic regime (c.f. Appendix D), 



where 



k ■ d r F g (r,k) = I g {r,k) 



I g (r,k) = - |n>(r,A0, F g (r,k) 
l q (r,p) = - I E>(r,p), F g (r,p) 



p-d r F q (r,p) = T q (r,p) , 



- |n<(r,fc), F g (r,k) + 1 

- j S < (r,p) , F q (r,p) - 1 



£122) 



(123) 



and the 'hatted' self-energy functions II, £ stand for 



8=1,2 



f^(r,p) = — \ u (P> s )^( r >P) u {P: s ) + v(p,s)Y^(r,-p)v(p,s) . (124) 

11 8=1,2 

The collision terms I g and I q on the right hand side of ( |122j ) describe the balance of 
gain and loss of partons in a phase-space element d s rd 4 k, or d s d 4 p, within a time slice 
around tq. Their explicit form is obtained as explained in Appendix D, and emerges as the 
result of applying the usual cutting rules [O] to the self-energies (|116|) , (|117|) and averaging 
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(summing) over initial (final) spin and color degrees of freedom (see Fig. 11). The resulting 
structure of the collision terms is as follows: 

(loss) i \ ~-(gain) 



bed 



/di^p^ f p f d^ p 

x | F a (l)F b (2) \M(ab -» cd)\ 2 9(q 2 ± - f gq ) [1 ± F c (3)] [1 ± F d (4)} 

- [1±F (1)][1±F 6 (2)] \M(cd^ab)\ 2 9(ql-n 2 gq ) F c (3)F d (4) } , (125) 

Here the notation is Ff(i) = Ff(r,pi) for the distribution functions of the parton species 
/ = 9iQiQ with four-momenta pi = p,P2,P3,P4 at space-time point r = (r , f). The 
structure of the collision terms in conjunction with the equations ( |122j ) is such that the 
squared matrix elements for the various scattering processes 12 — ► 34 (depicted in Fig. 11, 
and explicitly given in Appendix D) are weighted by a distribution function Ft{{) for each 
of the partons coming into the vertex and a factor (1 ± Ff>(J)) for each of the outgoing 
ones. The + sign is for gluons so that (1 + F g ) results in a Bose enhancement, and the — 
sign refers to quarks and antiquarks with (1 — F q ) indicating Pauli blocking. This is a direct 
consequence of the quantum statistical difference between the gluon and quark propagators 
(|72"|) and ([74]). The factors C a b(C c d) in front account for the identical particle effect, if 
incoming (outgoing) partons are indistinguishable. 

Eqs. ( |122j ) are the essential result of this Section. These Boltzmann-type equations are 
the final form of the transport equations for the dressed partons with phase-space densities 
F. The equations have a a drift term on the left hand side, and a collision term on the right 
hand side, which balances the various processes by which a dressed parton may be gained or 
lost in a phase-space element d 3 rd 4 p around time r°±Ar/2 = r ±l/(2/x). They describe the 
dynamics of the multi-parton system on kinetic scales, due to statistical, binary collisions, 
in which dressed partons appear as quasi-particles with a dynamical substructure, which 
is described in terms of probabilities to find a parton as a state consisting of a number of 
bare gluons and quarks of virtualities k 2 > fj, 2 q . These the underlying quantum fluctuations 
are embodied in F = M (g> V through the spectral density, or parton structure function, 



V and are determined by the renormalization equations (110). A scattering between two 



dressed partons is therefore described as a 'hard scattering' determined by the probabilities 
of finding in each of them a hard fluctuation with k 2 ~ q\ > fi 2 q of the order of the 
momentum transfer that sets the probing scale. This is expressed by the collision term 
on the right hand side, in which the products of F's involve the convolution of spectral 
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densities V, weighted with the squared matrix elements. A graphical illustration of this is 
shown in Fig. 12. 



3.6 A Monte Carlo calculation scheme 

A practical calculation scheme to compute the evolution of a multi-parton system as 
governed by the couled renormalization and transport equations may be outlined with the 
following concrete example. Consider the collision of two large nuclei with mass number 
A S> 1 at ultra-relativistic center-of-mass energy. Before their contact upon collision, 
the approaching nuclei appear as two highly Lorentz contracted discs of coherently bound 
gluons and quarks with a coherence length Lq = l///o, where fio ~ g s p± with p± = 
N g q/('KR\) sets the scale of the typical space-like parton virtuality, and hence Lq their 
characteristic transverse size. McLerran and Venugopalan have shown that if A is 
sufficiently large, the associated primeval parton distribution of the nuclei before, and 
shortly after the collision can be calculated non-perturbatively from first principles in terms 
of coherent quantum fields. Their conclusion is that, as long as the very early generation 
of this parton matter distribution has ApAr ~ 1, it cannot be described by a kinetic 
particle picture, which requires ApAr 3> 1. However, after a time to ~ l/z^o passed the 
nuclear contact, the parton matter has gone through a decoherence stage, so that the latter 
condition is be satisfied, and a kinetic description can be matched to the complex coherent 
evolution of the primeval matter. In other words, at time to, one may proceed with a 
probabilistic description of the parton dynamics in terms of the interplay between coherent 
radiative evolution and incoherent binary interactions, as suggested in the present work. 

At time to, one starts from the initial multi-parton state, and the subsequent time 
evolution of the partons' phase-space densities F gq {t, r,p 2 ,p) may be calculated by a Monte 
Carlo procedure, using the advocated discretization of space-time with four-dimensional 
cells of size Ar = AtA 3 r ~ /i~ 4 . 

(i) The first step consists in evaluating, from the collision kernel of the transport equations 
(122), the probabilities of scatterings among the initial partons within the time slice 



Ato = p 1 (to) = Mo 1 between to and to + Ato for each cell centered around r. The 
essential condition (112) provides the possibility of treating the scattering among 



partons incoherently, and requires that the impact parameter b a b of any two scattering 
candidates a and b must satisfy b^ < Pq 1 , implying for the momentum transfer of 
the scattering q^ b± > /j,q, where Pq is the initial virtuality of the partons at to, set 
by the coherence length of the colliding nuclei. The primary parton scatterings a + b 
that occur within Ato subject to this condition change the phase-space occupation of 
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partons at t± = to + Aio in two ways: on the one hand their spatial profile is altered 
due to gain and loss of deflected scatterers a' and b' in the spatial cells, and on the 
other hand, the virtualities are reset from k 2 = k 2 = $ to k' 2 = k' b 2 = q 2 bl _ > 

(ii) In the next step one must now calculate the quantum fluctuations, virtual plus real 
emission and absorption processes, that are triggered by the primary scatterings and 
the change of virtualities from k 2 to k 2 . That is, the parton structure functions 
need to be evolved within Aio according to the renormalization equations ( |110|) , 
which account for the associated renormalization and dissipation. One then obtains 
a spatial profile of new dressed partons, that defines the initial condition for the 
further evolution, starting at t\. The procedure repeats, as before at to, by evaluating 
the scattering probabilities in the next time slice At\ = between t\ and 

t\ + Aii, now however subject to the modified incoherence condition, that for any 
scattering of partons a and b their impact parameter must be b a b < max(A^, k 2 )^ 1 / 2 , 

It is important to realize that the condition ( |112| ) of a well defined separation between 
quantum and kinetic scales, imposes the crucial incoherence requirement for binary scat- 
terings, and allows a 'hard scattering' picture, in which the quantum evolution and hard 
scattering of evolving quanta factorize. The condition defines the range of validity for a 
probabilistic description in terms of incoherently scattering particles, and is essentially the 
uncertainty principle: the quasi-classical picture holds only, if the scattering partons may 
be treated as well distinguishable, incoherent quasi-particles of size l/y/k?, meaning that at 
least a 'formation time' of t^ ~ 1/yfk? must have been passed since their previous scatter- 
ing, during which their quantum structure evolves with virtuality k 2 . A concrete example 
of this scheme will be presented elsewhere. 

4. CONCLUDING REMARKS 

In this paper I have attempted to lay out a foundation to obtain from the fundamental 
quantum field theoretical principles of QCD a self-consistent kinetic description for the 
evolution of a high-energy system of self- and mutual interacting gluons and quarks. The 
main result is a set of two distinct, but coupled equations that govern the time evolution of 
the gluon and quark Wigner functions, the quantum analogues of the classical phase-space 
densities: 

(i) A renormalization equation, which describes the the momentum dependence of short- 
distance quantum fluctuations due to the partons' self- interactions. It defines the 
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state of a dressed parton as a quasi-particle with a renormalized mass and a decay 
width, corresponding to virtual and real emission and absorption processes. The 
solution of this equation describes, locally in space-time, the spectral density in terms 
of bare partons that are associated with the quantum substructure of the state of 
a dressed parton, and determines the partons' structure functions as well as their 
dynamical mass spectrum. 

(ii) A transport equation, which describes the space-time evolution of the dressed par- 
tons in the kinetic quasi-particle regime by means of mutual binary collisions. It 
determines the time dependence of both the change of the spatial density and the 
energy-momentum distribution of dressed partons due to elastic and inelastic colli- 
sions. Accordingly, it redistributes not only the partons in space, but also modifies 
their momentum spectrum and virtualities, which feeds back into the renormalization 
equations. 

Loosely speaking, the renormalization equation defines the state of dessed partons, whereas 
the transport equation governs the occupation of these states. The self-consistent solution of 
the equations provides the means to trace the dynamical development of the multi-parton 
system in real-time and full 7-dimensional phase-space d 3 rd 3 pdp 2 , spanned by position, 
momentum and invariant virtuality. It suggests a probabilistic, causal description, which 
is predestined for numerical evaluation by using Monte Carlo simulation techniques. 

The essential steps that lead to this kinetic framework may be summarized as follows: 

a) the path integral quantization of the classical action, using the CTP formalism with 

in-in boundary conditions, including initial state correlations at time to described by 
the density matrix /3(to); 

b) the 2-point source approximation, that allows to rewrite the CTP path integral as the 

generating functional for a possible color background field, and for the 2-point gluon 
and quark Green functions defined on the closed-time-path between to and too; 

c) the derivation of the self-consistent equations of motion for mean field (Ginzburg-Landau 

equation) and Green functions (Dyson-Schwinger equations); 

d) the transition to kinetic theory by imposing the physics-motivated well defined sepa- 

ration between the quantum scale that specifies the range of short-distance quantum 
fluctuations and the kinetic scale that characterizes the range of inter-particle corre- 
lations and stochastic binary interactions; 
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e) the conversion of the Dyson-Schwinger equations for the Green-functions to the set of 

renormalization and transport equations for the corresponding Wigner functions, on 
the basis of the separation of quantum and kinetic scales and a cellular space-time 
picture. 

f ) the calculation of the spectral density of bare partons, locally within the cells, from the 

renormalization equations, which defines the state of dressed partons in terms of their 
substructure, and the evaluation of the collision kernel of the transport equations, 
which determine the statistical occurrence of scattering processes among these dressed 
partons. 

This quantum kinetic framework may be extended in straightforward manner to include 
effects of a color background field, or gluonic mean field, that acts as a classical background 
medium in which the partonic quanta evolve (e.g. in a QCD plasma, where it may be gen- 
erated due to the bulk dynamics of soft gluon modes). This option has not been considered 
in the present paper, however the framework incorporates this possibility by considering a 
non- vanishing instead of setting it to zero as in Sec. 3. The inclusion of such a mean field 
would extend the set of renormalization and transport equations for the partons' Wigner 
functions, by coupling to the Ginzburg-Landau equation for the mean field. Qualitatively, 
the effect would be twofold: first, the poles of the Wigner functions would be shifted by 
a mean field generated mass /x(A^), and second, the transport equations would acquire an 
additional force term of Vlasov form. 

The future extensions and applications are manifold. Most important at first, I believe, 
is a detailed calculation for a specific situation where the concepts and formalism presented 
here may be illustrated and checked for consistency. For instance, it would be desirable to 
study a thermal (or close to thermal) parton system in this real-time description, and com- 
pare it with the well known results in the imaginary-time formalism of finite-tempaerature 
QCD. Such a project is planned to be carried out in the near future. On the other hand, the 
probabilistic interpretation of the real-time evolution of rather general multi-parton system 
offers the opportunity to simulate the dynamical development on the basis of the master 
equations with Monte Carlo techniques on a computer ]42| ]. 
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APPENDIX A: 

The CTP formalism and the 2-point source generating functional 



In this Appendix, I review the CTP functional formalism applied to the case of QCD. 
For additional reading concerning the general techniques, I refer to the extensive literature 
||, f|, H ||, [?], ||. In the in-in formulation of quantum field theory, mentioned in Sec. 2, 
the generating functional is defined as the in- vacuum to in- vacuum amplitude Z [J, p] = 
Tr X^(0 m |<^) j(ip\p\O m ) j, including possible initial state correlations represented by the 
density matrix p at to, and a sum over a complete set of states ip at too. hi the Heisenberg 
picture it is represented by 

Z P [J,p] = Tr^ P ewL^J p d*xJ f (xWf(x)\ p\ , (126) 

where / = g,u,u,d,d, . . . specifies the gluon and quark field degrees of freedom, and 4>f = 
(A^,^,tp). The symbol P refers to the time integration along a closed-time path in the 
complex t-plane as illustrated in Fig. 2: the path goes forward from to to too on the positive 
branch, and then back from too to to on the negative branch. The generalized time-ordering 
Tp is therefore defined such that any point on the second branch is understood at a later 
instant than any point on the first branch. Utilizing (|13|), eq. ( |126j ) can be rewritten as 

Z P [J + ,J_,p] = Tr{u\_(t ,t)Uj + (t,t )p(t )} , (127) 

where J+ (J-) is the source along the positive (negative) branch of Fig. 2a. In general 
J+ ^ J_, so that Zp depends on two different sources. If these are set equal, one has 
Zp(J, J, p) = Tr/i, which is equal to unity in the absence of initial correlations, being a 
statement of unitarity. The derivatives of Zp with respect to the sources generate the 
n-point CTP Green functions 

5JL\{x l )...5JL n n {x n ) } 

(128) 

where cti = ±, and the indices fi label the type of the i-th field. The functional Zp can 
be represented as a path integral by employing the relation between the Heisenberg and 
interaction pictures (|l6|). One imposes boundary conditions in terms of complete sets of 
eigenstates of the Heisenberg fields <&h at initial time t = to, 

<f> H (to,x)\<f> + {to)) = $i(t ,x)\<f> + {to)) = <j) + (x)\(j) + (t )) 
$tf(to,x)|<T(to)> = 0~(W(to)>, (129) 
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and in the remote future at t = t r 



c 



<p(x) | <p(too) ) . 



(130) 



Then, making use of the completeness of the eigenstates, one obtains from ( |127| ) the fol- 
lowing functional integral representation for Zp, 



Z P [J + ,J^,p] 



V4>+V(\>-Vlp ( 4>- (t ) | u]_ (t , too) I y(too) > 

X (<p(too)\ U J+ (too , t ) I </> + (to) ) ( + (to) I P I </> + (t ) 



;i3i) 



At this point it is convenient to represent + and — by greek indices a, /3, 7, . . . and to 
introduce a 2x2 matrix a as a 'metric tensor', 



.a/3 



a, /? 



t q /3 = = diag(l,-l) , 
and similarly higher rank tensors 

u Q/ 3 7 = 6 a pTpy , v af 3^s = sign (a) SapSpjS-fg , 



(132) 



(133) 



with the usual summation convention over repeated greek indices a, (3, ... With this con- 
vention one can generalize the classical action (Q) to account for all four field orderings on 
the closed-time path P: 



I 



I* 



+ l^[u aP ^ a A^\ u^d^A^AZ, v a ^ s A^AlAi) ,(134) 



where the correspondance with the terms I with the ones of @) is obvious (the color indices 
are suppressed here). Also, the following shorthand notation for the integration over the 
space-time variables will be used in the functional sense, 



J 4> = Id xJ(x) 4>(x) , 
Jp 



>Kd) = 



d 4 x<j)(x) K(x,y) <j)(y) . 



(135) 



Returning to the functional integral (131), I now utilize the above conventions and exploit 
the fact that the first two amplitudes are just the ordinary transition matrix elements in the 
presence of J+ and J_, whereas the density matrix element incorporates the initial state 
correlations at to at the endpoints of the closed-time path P, one obtains the path integral 
representation for Zp in analogy to usual field theory [|l0L O] 



Z P [J + ,J_,p] = J ]JV^ f exp i^I[</>?] + £ J /,^/j 



(136) 



Here DcpJ = DcpfDcj), , and I suppressed the formal presence of the Fadeev-Popov deter- 
minant associated with the gauge freedom, because for the class of ghost-free gauges (0) it 
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is equal to unity. The functional M[p] is the density matrix element containing the initial 
state correlations that may be represented by a non-local source functional K as follows: 

Af[p] = (<M* o )|p|0_(io)} = expLj^Kf^+A-]] (137) 
When expanded functionally as 

K f [4> + A-\ = K* + K&* + \k%$<$, + -K^f^fo}, + ... , (138) 



eq. (|136| ) becomes a functional of an infinite number of non-local sources M , which however 
contribute only at t = to, corresponding to the initial state correlations, 

Z P [J + ,J^,p] = Zp[J a ,K a/3 ,...] (139) 

where the constant term K* has been absorbed into the normalization and the local initial 
state kernel K* a {£) has been combined with the external source term jf(x), i.e. jf = 
4 + Kf. 

The corresponding generating functional for the connected Green functions Wp as given 
as usual by the logarithm of Zp, 

W P [J a ,K a P,K a ^,...] = -i lnZp[J a ,K af3 ,K a ^,...\ . (140) 

The functional derivatives of Wp with respect to the local sources Jf(x) define gluon 
and quark mean fields 4>J(x) as the expectation values of the single field operators, which 
can arise either through non-vanishing external sources, or, in the case of gluons, may be 
generated dynamically by the system itself depending on the initial conditions. Similarly, 
the variation of Wp with respect to the non-local kernels K^,(x, x'), K^?Zn(x, x' , x"), . . ., 
give the the n-point Green functions G^,(x,x'), G a ^7j,,{x,x' ,x"), etc., for the different 
particle species, which are the expectation values of products of n field operators. From 
(p9p and (p0|), one finds 



SW P _ ~ 

77T ~ ?f 

T77JT ~ 2 V ff ' 

0JX a/3 

8W P 1 



-gfp = ~ (Gf flfll + 3iGf f 4>}, + etc. (141) 

One immediately recognizes the inifinite hierarchy of the Green functions, the complete 
knowledge of which would correspond to the full solution of the theory. Clearly, in practice 
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one must truncate this infinite series. In what follows, I will assume that all n-point sources 
of order n > 3 (i.e. the kernels K a ^ etc., are neglegible and thus can be omitted. Such 
an approximation is justified if the higher order correlations generated by the n > 3 source 
terms are comparably small and the quantum dynamics of the system can be sufficiently 
well described by single-particle distributions corresponding to the 2-point functions. In 
fact, this is the very hypothesis of the parton description of QCD at large energies, where 
higher order correlations ('higher twist' effects) are kinematically suppressed by powers of 
a large momentum scale Q~ 2 corresponding to an approximate factorization of dominant 
short-distance 2-point correlations and larger distance 3-, 4-, ..-point correlations associated 
with multi-particle effects, an approximation which becomes exact in the asymptotic limit 
p6| . In this approximation the generating functional ( |139| ) then explicitly reads, 

iW P [J^,jJ,K^,k] 



Z P [J»,j,j,K^,k] = e 



T>A£D1>cPil> a exp i(l[A^ a M (142) 



which is the expression I stated in eq. ( pq ) of Sec. 2.2. Since the present interest concerns 
only cases where no external sources are present, one obtains from (|142| ) the mean fields 
and the 2-point functions for gluons and quarks by taking into account the fact in the 
that absence of external sources the establishment of a local colored mean field can only 
occur for the gluons, but not for quarks or antiquarks. Because of their bosonic character 
gluons the production of gluons can lead to a dynamically generated coherent field acting 
as a background medium, whereas the production of quarks and antiquarks occurs only in 
pairs and cannot yield a coherent mean field. Furthermore, a physical gluon mean field is 



determined the equality A^_ = A^_ = A^. Hence, one gets from ( 142 ), using ( 141 ) 



5W P ~ , . SW P SW P 



5J£(x) •■/'■■ 5j a (x) 6j Q (x) 



5W P 



SWp 



\ (iB#{x,y) + A^x)AUy)) 
iS a ?(x,y), (143) 



5k a/3 (x,y) 

where, as a reminder a, (3 = ±, not to be confused with color indices which I suppressed in 
this Appendix. 
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APPENDIX B: 
From quantum field description to kinetic theory 



In this Appendix I explain steps that lead from the Dyson-Schwinger equations (f42|), 
( f43| ) to the kinetic counterparts, the renormalization equations (^) and transport equations 
(|66|). A well separation between quantum scale of short-distance fluctuations and the 
statistical-kinetic scale, is the essential requirement for recasting the quantum theoretical 
problem, formulated in terms of the 2-point Green functions G(x, y), into the much simpler 
form of kinetic teory, employing Wigner transforms G(r,p). In the picture of cellular 
space-time, constructed in Sec. 3.1, a clearly defined separation between the two scales 
is controlled by the characteristic size of the space-time cells with volume Q ~ fj,~ , by 
choosing such that Ar qua < fx^ 1 < Ar^, i-e. the cell size is larger than the range of 
short-distance quantum fluctuations, Ar qua < /x~ 1 , but small compared to the mean- free 
path of the quanta between their kinetic, statistical interactions, Ar^m- This is illustrated 
in Figs. 6 and 7. Therefore the correlation between different cells is negligible by design, 
and the only relevant case is when the points x and y in the argument of G(x, y) lie in 
the same cell. In the interior of each cell, one can then assume approximate translation 
invariance, because large-distance inhomogenities of space-time are, by construction, not 
resolvable within the small cell volume. Thus, for each individual cell, one can Fourier- 
transform the Green functions over the cell volume, and use the common machinery of 
propagator theory as for homogenous systems, or the vacuum. Specifically, one transforms 
the Green functions G = D^, S with respect to the relative coordinate s = x — y, whereas 
the absolute coordinate r = ^(x + y) serves as a cell label: 

= / G (^'H = /!r ,! " IG( "> • (144) 

where G(r,p) is called the Wigner transform of G(x,y), and similarly for the self-energies 
£{x,y) = n MZ ,,£, 

= lw? e " q{ "" )£ (^ x - y ) - lw? e " , " £{r - a) ■ (145) 

If the separation between quantum and kinetic scales would be perfect (as in vacuum, 
where Arkin = oo), then the r-dependence would drop out and the Wigner transforms 
reduce to the ordinary Fourier transforms G(p), £{q). In case of moderately inhomogenous 
media, meaning that the Green functions and self-energies vary only slowly with r and 
are strongly peaked around s = x — y, as I consider here, then one can expand the Green 
functions G(x, y) = W (r + |s, r — = W(r, s) in a series of gradients, 

W(r + s,s) ~ W(r,s) + s ■ d r W{r,s) + ... , (146) 
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and analogous for the self-energies £. The basic Dyson-Schwinger equations 
remain the same in terms of the Wigner transforms, for example for the terms EG = Tl^D® 
or £G = US, one has 



d x'£(x, x') G(x', y) 



d 4 q d 4 p f 4 , 



(2tt) 4 (2tt) 4 



x + x' 



(147) 



-, x — X 



' j e -ip-(x'-y) Q 



x' + y , 
— n — , x -y 



The integrand will be significantly different from zero only if x' lies within the same cell as 
x and y, in which case \{x + x') ~ \{x' + y) ~ r. Therefore (|147f ) reduces to 

J d 4 x'£(x,x')G(x',y) = J -0^ | e -*f(*-») £ (r,p) G(r,p) + A(r,p)} , (148) 

where A(r,p) embodies the corrections to the ideal separation of cells. In terms of the 
gradient expansion (146), the first order correction which is 0(h), is given by 

i d£(r,p) dG(r,p) i d£(r,p) dG(r,p) 



A(r,p) 



2 dp p 



2 dr^ 



dpf 



[(dpE) ■ (d r G) - (d r £) ■ (d p G)\ 



(149) 



In general, the convolution between two functions / and g is given by 



j d 4 x'f(x,x')g(x',y) 



exp 



(<?(/). _ 9 (f). d ( 9 ) 



f(r,p)g(r,p) . (150) 



The gradient expansion ( |l4q) corresponds to keeping only the first two terms in the Taylor 
series of the exponential function, which gives the following set of conversion rules: 



d x'f(x,x')g(x',y) 
h(x)g(x,y) 
h(y)g(x,y) 
%f(x,y) 
d£f(x,y) 



f(r,p)g(r,p) + - [(d p f) ■ (d r g) - (d r f) ■ (d p g)) 
h(r)g(r,p) - - (d r h) ■ (d p g) 
h{r)g(r,p) + - (d r h) ■ (d p g) 
Hp" + \d?)f(r,p) 

+ ±d?)f(r,p). (151) 



If one applies these rules now to the Dyson-Schwinger equations ( fi2j) and (43), which upon 
setting for simplicity the mean field contributions fig = p q = 0, read, 

0^ D%{x,y) = 5 ab g^5 4 P (x,y) - d 4 x' 11^ ab , (x, x) Dy u h (x', y) 

D%(x,y) h y^P = S ab g^5 4 P (x,y) - J p d 4 x' L>£ aV {x, x') U^ b (x', y) (152) 
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and 



iT dx Sij(x,y) = Sij5p(x,y) + / d A x' T, ik (x, x') S kj (x', y) 



-Sij(x,y) ij- d y = Sij5 P (x,y) + / d x S ik (x, x) E kj (x', y) 



(153) 

one finds on the basis of the gradient expansion ( |146| ) a set of corresponding matrix equa- 
tions for the Wigner transforms of the gluon and quark Green functions, 

k 2 + \f T -ik-dr) D%(r,k) = 



d^(k)S ab i P - (I1D) 



ah 



1 ^2 



d^(k)5 ab 1 P - DU 



ab 



(d k U) ■ (d r D) - (drUn ■ {d k D»\ 



ab 

(154) 



(d k D) ■ (d r u) - (d r rf • (d k u] 



[MP 



ab 



and 



7 • (p + - d r ) ) Sij(r,p) = Sij 1 P + ( £ S ) + - ( (dpE) ■ (8 r S) - (d r E) • (d p S) 



Sij(r,p) 7 • (p 



Or 



5ij lp + ( S £ 



+ 



(d p S) ■ (d r E) - (a r 5) • ) (1,55) 



1/ 



(156) 



where 

1 forF, F 
for >, < 

recalling that G = D^ u , 5 and the self-energies £ = U^ u , £ each represent a 2x2 matrix 
as defined by (|27|). 



G 



G F G> 
G< G T 



£ 



£ F £> 
£< £ T 



(157) 



Adding the two equations of ( 154 ), respectively of ( [L55|) , yield the imaginary parts as the 
the passage to the renormalization equations stated in ([65]), whereas subtracting the two 



equations of (|154[) , respectively of (|155| ), gives the real parts as the transport equations 
(|66"|). For the gluon Wigner functions one obtains 



k 



1%) D£{r,k) 
k-d r D%(r,k) 



■dr(k)6ablp + - ( {U,D}_ 



1 r.v u (-) 



+ 7 ^ab 
ab 4 



fJIV 



7 ( [n,^L) + 



ab 4 

where [A, B\- = AB - BA, {A, B} + = AB + BA, and 



(158) 



g^i-) = \d^,d r x D° 



(159) 
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For the quark Wigner functions the corresponding equations read 



1 



- {7 -p, S ij {r,p)} + 

- {-f-d r , Sij(r,p)}_ 



5n 1 



ij i-p 



1 



- h-d r ,S]_ + - {E,S} 



IJ 



+ 1*P 
4 % j 



where 



2 ( h-p,s]_ 



2 ( 



(+) 



4 y 



(160) 



^+) = {«9 p A £,^s} -{d^,cf x s} . (161) 



The equations ( |160D for quark propagators can be formally brought in the same form as 
( |158| ) for the gluon propagators by multiplying the first equation of ( [L55| ) with 7 • (p + | d r 
)5u + Y^u from the left, and the second equation of (|155| ) by 7 • (p — | d r )5ji + E^; from the 
right, and then adding and subtracting the resulting equations: 



P 2 ~ \%) S, 



(7-p + s)^i P + i ( { s2 ' 5 } + ).. + I- 4 



4 'Sj 



8 V 



2 (7 • $o <y»j ip - I 



p • <9 r Sij(r,p) 
where the notation Ejj = <5jjE is employed, and 



S 2 , 5 



„4 (±) = I^( 7 . p + s) (^<-)+^(+)) ± (^(-)_^(+)) (7-p + s; 



B <±) 



1 



(7- <9 r ) (jrH + ± (JF(-) - ( 7 . d r 



(163) 



Due to the 2x2 matrix character of the equations ( |158| ) and ( |160| ), ( |162 ), the four compo- 
nents F,F,>,< of the Green functions D^ U ,S and self-energies IP^E mix, so that each 
of these equations actually represent a non-trivial coupled set of four equations. However, 



as advertised in Sec. 3.3, in the physical representation (73), 



G 



G A 
G R G c 



£ 



£ c £ R 
£ A 



(164) 



one has the great advantage that the retarded and advanced functions G R<yA ^ are deter- 
mined exclusively by the R and A components, and components, and only the equation 
for G c involves a mixing with these. Omitting for lucidity the gradient terms Q,F, the 
equations ( |158|) become in the physical representation a self-contained set for the retarded 
and advanced functions, 



k 2 



-di d 



'ab 



r, k) 



d^(k)6 ab + ~ ( U r Wd r W + D R WlL R W] 



\1V 

/ ab 



k-d r DZ R(A) (r,k) 



L f U R(A) D R(A) _ D R(A) n R(A)^ 



J ab 



(165) 
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and 



\ {tp, s$ A \r tP )} + = Sij - i ( [ 7 -a r , s»w]_) + \ ( {^ A \s^} 



I.I 



{7-d r ,sf/ A \r,p)} 



7 -p,S R ^ 



E R(A) ^ S R(A) 



(166) 



ij 



plus a set of mixed equations for the correlation functions, 



-d 2 \ rf u * ) < 



r, k) 



k-d r D%*(r,k) 



with D° v = D> u + £><„, and, 



[7-p, Sj(r,p)} + = 
i {-f-d r , S§(r,p)j + ■- 
with S c = S> + S<. 



1 

2 

2 











aft 










' ab 



j-p, S* 



+ (e*S^ + E*5^ - - S K £^(1Q8) 

jj V / ij 



APPENDIX C: 

The renormalization functions and the spectral densities of partons 



In the following I explain in more detail the steps that lead from the determining 
equation for the retarded self-energies IT and S, eqs. (p9|), (|100|), to the solution for the 



renormalization functions A g , A g , eqs. (104), and finally to the evolution equations for the 
phase-space densities F g , F q , ( 110 ). I exemplify the procedure for the simpler case of the 
quark self-energy. The case of the gluon self-energy is more eleborate, but conceptually 
it is completly analogous. I confine myself here to the leading log approximation (LLA), 
27] for additional reading. 

6ijp 2 t, eq. (H), given by 



referring to Refs. 



The quantity of interest is hence the quark self-energy Ejj 
(|95|). As explained in Sec. 2.1, when studying short-distance dynamics around the lightcone, 
it is appropriate and most convenient to work in the planar axial gauge n ■ A = 0, eqs. (§), 
(||) with the constant vector satisfying n 2 <C 1. Parametrizing it as r?/ 1 = (a+6, 0, 0, a— b) 
then requires n 2 = 4ab <C 1. Without loss of generality, one may set 6 = 1 and a <C 1/4, 
so that the scalar product of n with some four-vector q is n ■ q = q + + aq~ ~ q + with 
g ,± = go i 53, q + q~ = q 2 — q\ o± q 2 -C q +2 . Let me then proceed with eq. ( |100| ) for 
the variation of scalar quark self-energy function S(r, p) of a quark with momentum p and 
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virtuality p 2 <C p +2 , within a given space-time cell of volume Sl(r) = Ar°A 3 r = \i 4 (r) 
around r (c.f. Sec. 3.1): 

P 2 ^ R{A) (r, P ) = ~ 9l C F (2m) ** drT( P 'r) (169) 

Op 1 J (27T) 4 i J ro -l/(2fi) 

dA q (r,p') d 



x 



— 2 ( V 2 qg (r; p 2 ,p' 2 , k" 2 , k p , V , K k „) A g (r, fc")) X U™{j/, k" , n) 



dp' 2 dk" 

with 

Kp = ^=V +2 + \{P 2 -Pl) - P +2 (170) 
n z 2 

and 



w«v,a/» = -2 b / -n(7-n) + y-n ) ( 7 -n + 0(n2) (in) 

Here, p is the four-momentum of the incoming quark which branches into p' and k" of the 
outgoing quark and gluon, respectively. Because n 2 <C 1, the terms 0(n 2 ) are negligible 
and will be omitted in the following. 

The time integral J dr in (169) over the finite time slice = /i _1 (r) of the space-time 
cell around r, is weighted by the function T, satisfying 

/•CO 

/ drT(p'r) = 1, (172) 
Jo 

e.g. T = p'exp(— p'r), and the relation between momentum and space-time is determined 
by the uncertainty constraint, which limits the range of virtualities p 2 such that within 
the finite time slice Aro = fi^ 1 only those fluctuations k «-> k! + fc" are resolvable that are 
sufficiently short-living, with proper life-time tq ~ 1/k' and 7ro ~ k + /k'2 < Aro- (for 
details see Ref. Q) 

p+ 1 

t(p) = 7To(p) = — < t - ^ > (173) 
p z p(r) 

where To ~ 1/p' is the proper life-time of the virtual parton, to be understood in the 
averaged sense. 



To perform the integral (169), the procedure is as follows. First, due to the kinematic 
ordering condition p 2 ^> p 2 ,k 2 in the LLA, the gluon momentum k" can be decomposed 
as 



k" 



(l-z)p - (l-z) —p 2 + k'i + 0{n 2 ) 
V 2 / n ■ p 



?— k'i-p = k'l-n = . (174) 
p+ 



Then one can rewrite ( |171[) in the form 



Uf(p',k",n) = 2z{ 1 -p) + (l- z + 2z 2 )j^-p 2 . (175) 
y (n • k") 
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Next, one rewrites the the integration measure in (169) as 

' 2 k" 2 



d 4 p' = ^dp 2 dk" 2 dze\ ,r - '— 



.2 P 



z 1 — z 



(176) 



where the ^-function accounts for the aforementioned ordering of virtualities and acts as a 
kinematic constraint that limits the integration range, such that 

'2" 

Z 



k" 2 < (l-z)lp 2 - P — 



(177) 



Finally, to simplify the analysis, I approximate the r integral by [36] 



ro+l/(2jt) 



dr T(pr) 



1 



P 



(178) 



ro-l/(2/i) \P-(r) P 2 / 

where /i(r) is characterizes the size of the space-time intervall of the localized quantum 
fluctuations (c.f. Sec. 3.1, eq. fl48|). Using the above formulae, the integration over k" 2 
now gives 



xV* qg (r;p 2 ,A(l-z)(p 2 -r 



V (7 • n ) 
v (n • k") 

p +2 ,zp +2 ,(l - z)p +2 \ 



2z (7 -p) + (I- z + 2z 

2\ 



x A g r; (1 - z) \p z - 



2 P 2 \\ dA q (r,p') ( zp 



dp 



— ro - 1 



Next one integrates over p 2 , which yields for the last two factors in (|177| ) an 'effective 
vertex function' under the remaining z-integral 
2 

V 2 ff = ^V 2 g {r;p 2 ,zp 2 ,(l-z)p 2 ;p+ 2 ,zp +2 ,(l-z)p+ 2 ) A q (r; zp 2 ) A g (r; (1 - z)p 2 ) , 

(179) 



which in the LLA has been shown [26, 27] to generate the running of the coupling a s 
g 2 /(4ir), in the present case however modified by the finite time slize effect, 

n2 \ 



v; 



eff 



A-\r;p') a s ((1 - z)p l ) 9 I r - 1 



Employing this identification, inserting the decomposition 

2> 2 (7 • «) 



± R ( A \r,p) = Si( 7 -p) + S 2 



into ( |179| ) , and solving for A g and A in the parametrization of the self-energy (|9^) 



A,(r,p) 



1 



1 + Si + S 2 



A q (r,p) = A q (r;p 2 ) 



1 + Si 



(180) 



(181) 



(182) 
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one finds the form of the quark renormalization function as stated in Sec. 3.5: 

A q (r,p 2 ,p +2 ) = exp{- f ^ f X dzA{r,z,p 2 )^{z t e)\ , (183) 

where |^ 

At 2 \ a s{{ 
Alr.n .z) ~ — 

2tt \ p 



Afy.t,.) * (184) 



and 



with 



l + z 2 



^') = C F l _ . (185) 



2 2 2 

p n p z 



<P ) = 7T ^ = « 1 • ( 186 ) 

4(p-n) 2 p+ 2 

By repeating the same analysis for the gluon case one obtains 

A g (r,k 2 ,k +2 ) = exp|-£ + ^ j\zA{r,z^ 2 ) Q 7 f(*,e) + 7f(*,«0)J (187) 
where 

T?M = 2C * ( l- z + e(k+i) + 1 ^T + z{1 - z) ) ' 7 ^' £) = 5 ( z2 + ( 1 -^ 2 ) • 

(188) 

The renormalization functions Ay (/ = g,q) determine the form of the gluon and quark 
structure functions vj (r,p) defined in Sec. 3.1, eq. fl59|), as the spectral densities of the 
parton phase-space-distributions Ff(r,p): 

V f f '(r;x,p 2 ) = 6 f f ' 8(1 - x) S (f - /4) A f (r;n 2 gq ,p+ 2 ) 

+ A f {r,f gq ,p+ 2 ) Y;£til£dz [A(r, P 2 ,z) 7 ff\z,e)V f r (^,p 2 ;p +2 ) 

xAfL^p+^l . (189) 

Because Ff = Mf <8> Vf, eq. (plf), the variation of Ff(r,p) with the parton momentum 
p (more precisely, with virtuality p 2 ) therefore reflects the parton's changing gluon-quark 
substructure as dictated by the renormalization functions A g t q \, also called the Sudakov 
formfactor of a gluon (quark). This connection between Af and Ff emerges as follows. 
Treating gluons and quarks on the same footing, the differentiation of Aj 1 with respect to 
p 2 , the incoming partons's virtuality, yields 

d f 
p 2 —Aj 1 (r;p 2 ,p +2 ) = -Aj 1 (r;p 2 ,p +2 ) £ J dz A(r;p 2 , z) j f f 7 " (z,e) , (190) 



"'A detailed derivation of the effective coupling function A(r,p 2 , z) can be found in Ref. 
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where the sum over /" = g,q automatically fixes /' due to the symmetry properties of the 
kernels ^ under interchange of /' and f" |2(| . On account of the momentum sum rule 
fl82|) for the parton structure functions V(r,p) = V{r\ x,p 2 ), i.e. the momenta p + = xP + (r) 
of all partons in a given cell around r add up to the total cell momentum P + (r), one has 



1 = ^2 dx x Vj (r; x,p 2 ) , 
//' 

^ J dx x p 2 -^-^vj (r;x,p 2 ) = 
and therefore can rewrite (|185|) in the following form: 



and hence 



(191) 



(192) 



a 



V f f '(r;x,p 2 ) + A f 1 (r;p 2 ,p +2 ) 



f 



f 



1^. nr? ^+ 2N 



2 ^ T>f ( 2\ 

V -q^P) (r;x,p ) 



= -A]\r;p 2 ,p +2 )^2 J dz A(r-p 2 ,z) 1 f f ' f "(z,e)^V f r (r; ^,zp 2 )(193) 

Employing eq. ( |19C| ), then yields 

Af(r;p 2 ,p +2 )p 2 —,V f f '(r;x,p 2 ) = (194) 



f x 2 

■i„ [r;-,zp 
z 3 \ z 



A-\r;p 2 ,p +2 ) I dz A(r;p 2 , z) f '' \z,e) Ivf (r;x,p 2 ) - -V 
The final evolution equation for the parton phase-space densities Ff(r,p) is obtained by 



(i) multiplying with Af, (ii) convoluting the resulting equation according to eq. (56) with 
the local parton density Nf(r,p), i.e. the number of dressed partons in a given cell around 
r, and, (iii) accounting for the competition between real emission and reverse absorption 
processes pTj], using the fact that the squared matrix-elements oc are invariant under 
reversal a — > be and be — > a. The extended result is: 



where 



1 

-F g [r;x,k 



f 1 dz A{r;k 2 ,z) | 

+ 2N f F q (r;x,k 2 )T^(z,e) - N f F g (r; x, k 2 ) Yf (z, e) 



-F g (r-A,zk 2 
z V z 



rf(z,e) 



dz A(r;p 2 ,z) j ^ F q (r; ^, zp 2 ) - F q (r;x,p 2 ^) 
+ F g (r;x,p 2 )Tf(z,e\ 



ft 



Ff 
F f , ± 1 



(195) 



(196) 
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represents the net emission probability, being a manifestation of the principle of detailed 
balance between the rate of emission of parton f from a parton /, and the rate of absorption 

£l £ll 

of a quantum / in the phase-space proximity of parton /. The net rate J results in a 

suppression when F g or F q becomes large, and thus reflect correctly the Bose-Einstein and 

Fermi-Dirac statistics (+ for gluons, — for quarks). On the other hand, when F g and F q 

f f" 

are small compared to 1, the usual branching kernels are recovered. 



APPENDIX D: 

Derivation of drift term and collison kernel of the transport equations 



This Appendix explains the derivation of the transport equations (|110| ), ( |122|) that 
govern the kinetic, dispersive dynamics of dressed gluons and quarks. First of all, one 



observes that from the matrix representation of the CTP Dyson-Schwinger equations (37), 
( |38|) in terms of the Green functions G F , G F , G > , G K and self-energies £ F , £ F , £> , £ < , where 
G = D^ U ,S and £ = 11^,2, follow immediately the corresponding equations for the 
retarded, advanced, and correlation functions G R ,G A ,G C and £ R ,£ A ,£ C , as given by 

(I 



D 



-i 



[MP 



R(A) J k ) 
S R(A)( r > P ^ 



D 



-l 

(O)R(A) 



n 



R{A) 



and 



Sc{r, p) 



D R 



S R 



S (0)R{A) 



D {<S)C 



J R{A) 



(197) 



D (0)C 



S 



c 



S A 



(198) 



for the quarks. It is convenient to introduce scalar and dimensionless self-energy functions 
II and E through 



n£(r,fc) = Sab (k»k? - g^k 2 ) U(r,k) 

so that the propagators can be written as 

1 



Sy(r,p) = 5 llP 2 ±(r,p) , (199) 



D% RiA \r,k) = 5 ab (-d, v {k)) 



S* (A \r,p) = Sij (7 -p) 



1 



cr ± (Ti ' 



(200) 



and the correlation functions as 



D 



C 

ab 



(r,k) = S ab (-d^(k)) 



n£ + vrf 



[1 + F g ] &(r,p) = Sijfr-p) 



-2ia 2 



(T ± CTl 



[1 " F q ] 



(201) 
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where d^ u (k) is defined in (||), and F g and F q are the phase-space densities of dressed glu- 
ons and quarks as defined by (|56|), This separation of real and imaginary contributions 
uniquely determines the Green functions in terms of the three functions 7Tj, respectively o"j, 
as can be shown rigorously Q. The real parts correspond to the dispersive and wavefunc- 
tion renormalization piece, whereas the imaginary parts give rise to dissipation and decay. 
Formally, 

vr = k 2 (1 - Reft) Reft = - (lL R + fl A 
7f : -k 2 Imtl Imft = - (fl R - IT 4 ) (202) 



i 
2 

a = p 2 (1 - ReS) ReS = ^ (t R + £ y 



rr, - -p 2 ImS ImS = - (£ R - £ A ) , (203) 



and 



7r 2 = | (ft< + fi>) a 2 = 1 (S< + £>) . (204) 

Next, recalling that the correlations among different dressed partons determine their mutual 
interactions at kinetic space-time scales, one focuses on the correlation functions and 
S c . Noting that D A = D R J , S A = S R ^ , and employing the representations flSSj), 

D°,(r,k) = -2m (-d^(k)) [1 + 2F g (r,k)] 8 (k 2 - M 2 g (r,k)) 
S c (r,p) = -2ni ( 7 -p) [1 - 2F,(r,p)] 8 (p 2 - M 2 q (r,p)) , (205) 
where color indices are suppressed, eqs. fl200|) - (|204| ) may be combined to write 



-cT(k) D%{r,h) = (vr + ii^y 1 H g - H g {it Q - ivn)" 1 = -2% 2 ^ 2 

vr + 711 

~1-pS c {r,p) = (cjo + ivir 1 H q - H q {a - iai)~ l = -2i 2 ^ (206) 

°o + CT i 

where (c.f eq. (85)) 

H g (r,k) = 1 + 2F g (r,k) H q {r,p) = 1 - 2F q {r,p) . (207) 

Hence, on account of ( |197| ) and ( |198| ) one finds 

-d^(fc) (Z)^- 1 - n C )^ = -(i/ g VT - TT H g ) + i {H g TT\ + 7T1 # s ) 

7-P (^gf 1 " SC ) = -(#9 ff o - #,) + i (H q <n + (7i #,) . (208) 
Then, by inserting the expressions Q203;), (p04j) for tti and cr,, one obtains 



F g n -7r F g = i((i + i^)n< + n< (1+F g ) - F g fl> - n>F g ) = c 9 
F 9 a -a F g = -if(l-F ? )S< + S< (1 - F,) + F g S> + £> F q ) = C q , (209) 
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with the self-energies II < and S < given by ( |116| ) and ( |117| ) together with ( |199| ). Finally, as 
argued in Sec. 3.5, the presumed clear separation of quantum and kinetic space-time scales 
allows one to treat the kinetic dynamics quasi-classically, by expanding both sides of ( |209| ) 
in terms of fi and keeping only the lowest order contributions. The lowest order non-zero 



terms on the right hand sides of (209) correspond then to the Born collision terms which 
are of order h, 

C g (r,k) = (1 + F g ) n< - F g IL> + O (h 3 ') 

C q (r,p) = -(1 - F q ) t< + F q t> + 0(h 3 ) , (210) 



whereas the lowest order non-zero terms on the left hand sides of ( |20S| ) result in the classical 
Poisson brackets which are also of order h, 

F g ttq — ttq Fg^j (r, k) = -i ( (d r F g ) ■ (d k ir ) - (d k F g ) ■ (d r n )\ 

k ■ d r F g (r, k) - \d r k 2 Q ■ 8 k F g (r, k)) + 0(fi 3 ) 



( dnojr, k) 
V dk 2 



•j 'i 



(^F q a - a F q ^j (r,p) = -* ( (d r F q ) ■ (d p a ) - (d p F q ) ■ (d r a fj (211) 
'da (r,p)' 



dp 2 



p-d r F q (r,p) - -d r p 2 -d p F q (r,p)) + 0(fi 6 ) 



where as before d r = d/dr^, d k = d/dk^, etc., and the dot denotes a scalar product 
of four-vectors. The latter equalities in these two equations are obtained by using the 



fact that the solutions of the dressed partons' energy spectra ( 119 ) are strongly peaked 



around momentum transfers q\ ~ fi 2 q , i.e. ko(r, k) ~ y k 2 + fi 2 q and po(r,p) ~ yp 2 + n 2 q , 
because of the well known QCD specific logarithmic behaviour of the spectral densities 
oc a s \n(qj_/ij,g q ), and the power law form of the scattering cross-sections oc a 2 qj n (n ~ 4). 

Finally, using d r ko = d r po « 0, and equating ( |209D and ( |2 11 ) , one obtains the transport 
equations of Boltzmann type, stated in Sec. 3.5, eqs. (|115|), 



k-d r F g (r,k) = l g (r,k) p-d r F q (r,p) = l q (r,p) , (212) 

where the Lorentz invariant collision terms I on the right hand side are defined by 



V dk 2 

( dao(r ,p) 
\ dp- 



k 2 =fi? gg 



C ^P) = -M^^) Z ^P)- (213) 

p 2 =v 2 9q 



In (pll[) and (|213|), the derivatives with respect to the virtuality k 2 are to be taken at [i 2 , 



~gqi 

which, according to (|54[) and (p5), defines the scale at which a dressed parton appears as 
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a quasi-particle, with the renormalization and dissipation effects taken into account in the 
spectral densities V g and V q , eq. ( |109|) . Their form is determined by the renormalization 
equations ( |110| ), and therefore the parton distributions F = M ® V contain implicitely the 
short-distance quantum effects. The derivatives of 7To and cxo are related to the renormal- 
ization functions A g and A q , respectively, via the correspondence of the representations 
( |200| ) and (p^), and one finds 



V dk 2 

-(^L_., = W^)\^„ = 1, (214) 



where the latter equality results from the normalization condition (|94|), 

The explicit forms of the collision integrals X is obtained by substituting the correlation 



functions ( |205| ) into the two-loop expressions (|116| ), ( 117 ) for the self-energies II, X, and then 
inserting those into eqs. ( |209| ) . Applying the standard cutting rules [^] to the resulting 
self-energies, as symbolically represented in Fig. 11, yields the different binary collision 
processes ab <-> cd by which a parton of type a may be gained or lost in a phase-space 
element, namely gg <-> gg, gg «-> qq, gq <-> gq, qq <-> qq, qq <-» qq. The corresponding 



collision integrals Z a may be compactly represented in the generic form of (|12 



bed 



n f d 3 p 2 f rf 3 p 3 f d 3 p 4 4 4 

= -JLCabC cd j J^y^J J^^Jj^El^) 6 (P1+P2-P3-P*) 

x | F a (l)F b (2) \M(ab - cd)| 2 %i - $ q ) [1 ± F c (3)] [1 ± F d (4)] 

- [l±F a (l)][l±F 6 (2)] |A4(cd^a6)| 2 %i-/i^) F c (3)F d (4) } , (215) 

where the F a (i) = F a (pi,r) denote the distribution functions of parton species a = a,b,c,d 
and corresponding four-momenta pi = Pi,P2,P3,P4 at space-time point r = (r°,r). As a 
consequence of the representations (205), the squared matrix elements \M\ 2 for the pro- 
cesses ab «-> c<i (which contain the 2 — ► 2 kinematics, color and spin structure, as given 
below) are weighted by a distribution function F a for each of the particles coming into 
the interaction vertex and a factor [1 ± F a ] for each of the outgoing ones, with the + sign 
refering gluons and the — sign to quarks and antiquarks. The factors S a b = (1 + ^ab)" 1 and 
S c d = (1 + cW)" 1 account for the cases where the two incoming and/or outgoing partons 
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are identical. Using the identities 11 , 

d^u(k) = ^2 Efj,(k,s) ■ e* v {k,s) 

s=l,2 

7-p = 2p u(p, s) u(p, s) = 2p ^2 v(p,s)v(p,s) , 



(216) 



5 = 1,2 



s=l,2 



finally, the squared matrix elements are obtained by evaluating the amplitudes illustrated 
in Fig. 11, squaring those, averaging over initial colors and spins, summing over final colors 
and spins, and summing over quark flavors. The resulting expressions are standard and 
given by: 



\M(g a gb -> g c gd)\ 2 = 



9i 



(8 • 2)< 



E 



9tt' faedfebc \ptg( \ \t' ' uv i \ 

(p 1 _ p4 )2 A (-Pl)Pl -P4,P4j A ^ {P2-P3,-P2,P3) 



color, spin 

. 9tt' faecfebd \ryrvi \ \r' ua I \ 

+ ~ ~ A (-Pl,pi -P3,P3) A " (p2 -P4,-P2,P4J 



+ 



(Pi - P3) 2 
9tt' fdbefecd 
(Pi +P2) 2 



(217) 



A PMT (-pi, ~P2,Pl + P2) A T " CT (-p 3 ~Pa,P3,Pa) 

2 



+ V P a bcd(P^P^-P3,-Pi) 



\M(g a g b -> &<?j)| 2 



gs y y 



%(P4) ( ^7 • £ fe) 7 (p | P ! J 2 3) 7 • £(P1) ) 1*(ps 



(P2 - P3) 2 

+ «i(P4) ( i / g6C;Z 0- £ft ^ ) +ff ^ Vp(-P1.~P2,P1 +P2)) «i(P3) 



(218) 



(8-2)(3-2) 



E E 



/cabrpc 



(pi +P2) 2 



7 yT~ \uv(Pl ~P3, ~Pl,P3) Uj(p 4 h P Ui(P2) 

(Pi - P 3 ) 2 



- iT^Tgj % (p 4 ) 7 • £(pi) 7 / (P2 ^ 7 • e(P3) «i(?2) 

(P2 - P3j 

- iZ£Zg % (p 4 ) 7 • e»(p 3 f {P2+P l\ • £ ( P1 ) u,(p 2 ) 

IP2+P1J 



(219) 



\M(qiq k -> ^g/)| s 



.4 



(3 • 2)* 



E E 

/11/2 color, spin 



2-.r fe > j (p 4 )7^ i (p 1 ) 



(pi - p 4 ) 2 



MP3h^ u k(P2) 



11 Note that in contrast to the standard normalization for fermions oc yjm/p°, here the normalization is 
chosen commonly for both gluons and quarks oc l/(2p°). 
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- S flf2 T il T kj u j(P4,huU k (P2)-, v2«i(P3)7i/«i(Pl 

{Pi - P3) 



(220) 



\M(qiq k -> qjqi) 



9t 



3-2-2 



E E 



5 hu 6 hh T z a j T i a k u j(Piht^ u i(pi 



/11/2 color, spin 
~ 5 hh 5 Uh T ik T liUj{Pi)lvVl{p?,) 



Ol - P4) 2 



1 



(Pi -P3) 



Here 



(221) 
(222) 



A^(P1,P2,P3) = (P1"P 2 )V P + (P2-P3)V + (P3-Pl)V 
Ol>P2,P3,P4) = Zafce/cde^V 7 " S'V) + facefbde^g™ ~ lT <T) 

+ fade fcde (9^9^ ~ <f»<T) , (223) 

are the usual 3-gluon vertex function, and the 4-gluon vertex, respectively. The shorthand 
notation suppressing spinor and polarizarion indices, u(pi) = u(pi,si) a , e(p2) = z(P2,S2), 
etc., is employed, and in Q22C| ) and ( |22lD , 6ff is equal to 1, if the flavor of the two quarks 
are of the same flavor, and is zero otherwise. 
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FIGURE CAPTIONS 



Figure 1: 

Illustration of the difference of expectation values in the in-out and the in-in formalism, 
corresponding to the time-ordered product of field operators. 

a) In the usual S'-matrix formalism with a trivial (diagonal) density matrix p(to) = 1 
and |0 m ) = |0 OU *), it suffices to calculate (0 OU *| . . . |0 m ), because of the symmetry of the 
time paths (to, too) and (too, to). 

b) In the general case of a non-trivial initial state with multi-particle correlations de- 
scribed by p(to) 7^ 1, one must account for the complete time evolution on a closed-time-path 
from to to too and back to to by calculating (0 m | . . . p\0 m ). 

Figure 2: 

a) The close-time-path in the complex t-plane for the evolution of operator expectation 
values in an arbitrary initial state. Any point on the forward, positive branch to — > too is 
understood at an earlier instant than any point on the backward, negative branch too — ► to- 

b) The four different possible time orderings (ti^t^) in the arguments of the 2-point 
Green functions G(x,y) = G(t\,x;t2,y), corresponding to G F ,G > ,G < ,G F . 

Figure 3: 

Matrix representation of the CTP 2-point functions: a) The Green function G(x,y), 
and b) the self-energy function T,(x,y) oc [G(x,y)] 2 . 

Figure 4: 

Diagrammatic representation of the Green functions (i) G^(x,y), the of the bare prop- 
agators, (ii) Gro)(x, y), including the effect of a mean field by dressing the bare propagators 
with a dynamical mass, and (iii) G(x,y), the full propagators, dressed by both local mean 
field and non-local quantum self-interactions (Dyson-Schwinger equations). 

Figure 5: 

a) Diagram of the function Tp , eq. (j34|), representing the sum of all two-particle 
irreducible graphs of order h 2 , ffi , . . ., with fully dressed propagators D^ u and S. 

b) Illustration of the self-energies IP^ and S, eqs. ([!(]) and (41), which derive from 

(2) 

T p by functional differentiation with respect to and S. 
Figure 6: 

a) Classification of the different scales of relevance: (i) the quantum scale Ar qua , of 
the order of the spatial extent of quantum fluctuations associated with the 'radiative' self- 
energies, and defining a dressed parton state as a quasi-particle; (ii) the kinetic scale Ar^n, 
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measuring the range of correlations and binary interactions between these quasi-particles, 
with the 'collisional' self-energies; and (iii) the 'macroscopic' scale Ar mac , where the dy- 
namics can be described in terms of bulk thermodynamic variables or hydrodynamics. 

b) The quality of separation of quantum and kinetic scales is controlled by the choice 
of //(r). Because the the multi-particle dynamics of the system in general may change 
the scale of separation in space-time, one may choose fj,(r) variable to optimize the kinetic 
description. 

Figure 7: 

a) Illustration of the cellular space-time picture, with cell size chosen intermediate 
between quantum and kinetic scales such that the separation between the two scales is 
optimal, so that short-distance quantum correlation between different cells are negligible. 

b) Representation of the partons' phase-space densities F = N <8> V as a convolution 
of the statistical density of dressed partons N with the spectral density V of each dressed 
parton, describing its intrinsic density of bare parton states as its quantum substructure. 

Figure 8: 

In the cellular space-time picture, the 'absolute' coordinate r labels the kinetic space- 
time dependence 0(Ar/ t i n ), whereas the 'relative' coordinate s measures the quantum space- 
time distance 0(Ar qua ). 

Figure 9: 

The 'radiative' self-energies in one-loop approximation, eqs. (|95|): a) the retarded 
(advanced) gluon self-energies IL^j^ , b) the retarded (advanced) quark self-energies T, R ^ . 



Figure 10: 



The 'collisional' self-energies in two-loop approximation, eqs. (|116j) and (117): a) the 



contributions to the gluon correlation functions II< , b) the contributions to the quark 
correlation functions E^. 

Figure 11: 

Cutting the 'collisional' two-loop self-energies, gives the different binary 2 — ► 2 collision 
processes, namely a) the gluon terms gg <-> gg, gq «-> gq, gg <-> qq, and b) the quark terms, 
qg qg, qq <-+ qq, qq <-> qq. 
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Figure 12: 

Illustration of the 'hard scattering picture', for the evolution of a multi-parton system 
on the basis of the coupled renormalization and transport equations: 

a) A dressed parton is described as a quasi-particle with a dynamical substructure, 
corresponding to an instantanous state consisting of a number of bare gluons and quarks 
(its radiative cloud). These the underlying quantum fluctuations are embodied in the spec- 
tral densities, or parton structure functions, which are determined by the renormalization 
equations. 

b) A binary collision between two dressed partons is described as a statistically occurring 
'hard scattering', determined by the local density of dressed partons, and convoluted with 
their spectral densities at the 'hard scattering scale' of the order of the momentum transfer. 
This is described by the transport equations. 
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